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Abstract 



In this thesis I present a wide-ranging study of locahsing relational degrees of 
freedom, contributing to the wider debate on relationism in quantum mechanics. A 
set of analytical and numerical methods are developed and applied to a diverse range 
of physical systems. Chapter El looks at the interference of two optical modes with 
no prior phase correlation. Cases of initial mixed states — specifically Poissonian 
states and thermal states — are investigated in addition to the well-known case 
of initial Fock states. For the pure state case, and assuming an ideal setup, a 
"relational Schrodinger cat" state emerges localised at two values of the relative 
phase. Circumstances under which this type of state is destroyed are explained. 
When the apparatus is subject to instabilities, the states which emerge are sharply 
localised at one value. Such states are predicted to be long lived. It is shown that 
the localisation of the relative phase can be as good, and as rapid, for initially mixed 
states as for the pure state case. Chapter El extends the programme of the previous 
chapter discussing a variety of topics — the case of asymmetric initial states with 
intensity very much greater in one mode, the transitive properties of the localising 
process, some applications to quantum state engineering (in particular for creating 
large photon number states), and finally, a relational perspective on superselection 
rules. Chapter HI considers the spatial interference of independently prepared Bose- 
Einstein condensates, an area which has attracted much attention since the work 
of Javanainen and Yoo. The localisation of the relative atomic phase plays a key 
role here, and it is shown that the phase localises much faster than is intimated in 
earlier studies looking at the emergence of a well-defined pattern of interference. A 
novel analytical method is used, and the predicted localisation is compared with the 
output of a full numerical simulation. The chapter ends with a review of a related 
body of literature concerned with non-destructive measurement of relative atomic 
phases between condensates. ChapterElexplores localising relative positions between 
mirrors or particles scattering light, addressing recent work by Rau, Dunningham 
and Burnett. The analysis here retains the models of scattering introduced by those 
authors but makes different assumptions. Detailed results are presented for the 



case of free particles, initially in thermal states, scattering monochromatic light and 
thermal light. It is assumed that an observer registers whether or not an incident 
light packet has been scattered into a large angle, but lacks access to more detailed 
information. Under these conditions the localisation is found to be only partial, 
regardless of the number of observations, and at variance with the sharp localisation 
reported previously. 
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Introduction 



1.1 Preface 

It is a widely accepted principle of modem physics that absolute physical quantities 
have no intrinsic usefulness or physical relevance. Typically however, it is difficult 
to determine the extent to which some particular theory can be explicitly formu- 
lated in relational terms. There is an ongoing debate, and a substantial literature, 
which explores different aspects of quantum mechanics from a relational point of 
view, within diverse fields including quantum information, quantum gravity and 
foundational studies of quantum mechanics. This thesis contributes to this activity 
by presenting a comprehensive study of systems wherein some relationally defined 
degree of freedom "localises" — becoming well-defined, exhibiting strong correla- 
tion, and becoming in some sense "classical"^ — under the action of some simple, 
well-characterised dynamical process. Three physical systems are studied in depth, 
each built upon a simple measurement-based process. The emphasis is then on the 
induced, post-measurement properties of states of these systems. Chapters El and El 

^The precise meaning of "classical" here depends on the physical system being studied. Con- 
sider, for example, an optical system wherein a number of "quantum" light sources are phase-locked 
by processes of localisation, and are subsequently fed into a system of phase shifters, beam-splitters, 
and detectors. The dynamical properties of the light are predicted to be the same as for "classi- 
cal light" fields with corresponding intensities and phase differences. These fields are generically 
described by complex numbers and have absolute optical phases. 
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consider the localisation of relative optical phase in interference experiments wherein 
light from independent sources leaks onto a beam-splitter whose output ports are 
monitored by photodetectors. Chapter |3] looks at the spatial interference of inde- 
pendently prepared Bose-Einstein condensates, a process wherein the localisation 
of the relative atomic phase plays a key role. Finally, Chapter discusses the lo- 
calisation of relative positions between massive particles as they scatter light. This 
thesis contributes many new results using both analytical and numerical methods. 
Up to now, there has been little attempt to develop in detail the issues common 
to examples such as these. This thesis lays out a "modus operandi" that can be 
applied widely. Much of the content of this thesis has been published in |Cable05] . 

As an example of the relevance of this thesis topic, consider the highly involved 
controversy concerning the existence or otherwise of quantum coherence in several 
diverse contexts — some key examples being the quantum states of laser light, Bose- 
Einstein condensates and Bardeen-Cooper-Schrieffer superconductors. For a recent 
overview and fresh perspective on the controversy see |Bartlett?)5] . To take just one 
example, M0lmer in his well-known publications |M0lmer97al |M0lmer97b| challenges 
the assumption, common in quantum optics, that the state of the electromagnetic 
field generated by a laser is a Glauber coherent state, with a fixed absolute phase 
and coherence between definite photon numbers. By treating carefully the gain 
mechanism in a typical laser, and pointing to the absence of further mechanisms to 
generate optical coherence in standard experiments, he concludes that the correct 
form for the state of a laser field is in fact a Poissonian improper mixture^ of number 
states with no absolute phase and no optical coherence. A key question is then why 
two independent laser sources can demonstrate interference, as has been observed 
experimentally. M0lmer argued that in fact there is no contradiction, by detailed 
numerical studies. He showed that a suitable process of photodetection acting on 
optical modes which are initially in number states, can cause the optical modes to 
evolve to a highly entangled state, for which a stable pattern of interference may 
be observed. This corresponds to the evolution of a well defined correlation in the 
phase difference between the modes. Chapter E] takes up this example. 



^An improper mixture is a density operator which cannot be given an ignorance interpretation. 
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In studying localising relative degrees of freedom a number of questions suggest 
themselves regardless of the specific physical realisation that is being considered. 
How fast are the relative correlations created, and is there a limit to the degree lo- 
calisation that can be achieved? What are the most appropriate ways of quantifying 
the degree of localisation? How stable are the relative correlations once formed, for 
example against further applications of the localising process with additional sys- 
tems, interaction with a reservoir, and the free dynamics? Does the emergence of 
relative correlations require entanglement between the component systems? Does a 
localised relative quantum degree of freedom behave like a classical degree of free- 
dom, particularly in its transitive properties? What role does an observer play in 
the process of localisation? 

In addressing questions such as these, this thesis has several key objectives. 

In addition to revisiting the more commonly considered examples of pure initial 

states,'^ the focus is on initial states which can readily be prepared in the laboratory 

or are relevant to processes happening in nature, and in particular on examples of 

initial states which are mixed. Working with mixed states it is important to be 

wary of common conceptual errors and in particular of committing the preferred 

ensemble fallacy. The fallacy is to attribute special significance to a particular 

convex decomposition of a mixed state where it is wrong to do so. Another key 

goal is to simplify the analyses as much as possible. The emphasis is on deriving 

analytical results rather than relying exclusively on stochastic numerical simulations. 

A further goal is to identify descriptions of the various measurement processes as 

positive operator- valued measures (POVM's), so as to separate out the characteristic 

localisation of the relevant relative degree of freedom from the technical aspects of 

a particular physical system. Identifying the relevant POVM's can also facilitate 

analogy between localisation in different physical systems. And finally there is a 

preference for specifying in operational terms preparation procedures, and measures 

■^The studies of localising relative optical phase in |M0lmer9'7al |M0lmer97b| ISandersOS] con- 
sider specifically the case of initial photon number states, which are challenging to produce ex- 
perimentally (in fact M0lmer's numerical simulations assume initial Fock states with order 10^ 
or 10^ photons). Most of the existing literature concerning the spatial interference of inde- 
pendently prepared Bose-Einstein condensates assumes initial atom number states, for example 
j,Tavanainen96bl IYoo97j . The analysis of localising relative positions between initially delocalised 
mirrors or particles in |Rau03 .Dunningham04, assumes initial momentum eigenstates. 
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of the speed and extent of localisation, rather than relying on abstract definitions 
with no clear physical meaning. 

At this point mention should be made of some further topics with close connec- 
tions to this thesis. The first is that of quantum reference frames. Put simply, a 
frame of reference is a mechanism for breaking some symmetry. To be consistent, 
the entities which act as references should be treated using the same physical laws 
as the objects which the reference frame is used to describe. However in pursuing 
this course in quantum mechanics there are a number of immediate difficulties to 
address. There is the question of the extent to which classical objects and fields are 
acceptable in the analysis. Furthermore, in quantum mechanics establishing refer- 
ences between the objects being referenced and the elements of the reference frame 
causes unavoidable physical disturbances. Careful consideration of the dynamical 
couplings between them is necessary, and in particular the effects of back-action 
due to measurement. Finally, translation to the reference frame of another observer 
itself requires further correlations to be established by some dynamical process (in 
contrast to the simpler kinematical translations between frames possible in classical 
theories) . 

Another recurring topic is that of superselection rules. In the traditional ap- 
proach a superselection rule specifies that superpositions of the eigenstates of some 
conserved quantity cannot be prepared. This thesis involves several examples where 
modes originally prepared in states (pure or mixed) with a definite value of some 
conserved quantity express interference, but without violating the corresponding 
superselection rule globally. Following the more relational approach to superselec- 
tion rules of Aharonov and Susskind Aharonov67j states thus prepared, with a 
well-defined value for the relative variable canonically conjugate to the conserved 
quantity, may be used in an operational sense to prepare and observe superpositions 
that would traditionally be considered "forbidden". There is also the question as 
to whether superselection rules contribute to the robustness of the states with a 
well-defined relative correlation. In particular, if typical dynamical processes obey 
the relevant superselection rule then averaging over the "absolute" variables does 
not affect the longevity of these states. 
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1.2 Executive summary 

Chapter [5J Localising Relative Optical Phase 

Chapter 121 looks at the interference of two optical modes with no prior phase cor- 
relation. A simple setup is considered wherein light leaks out of two separate cavities 
onto a beam splitter whose output ports are monitored by photocounters. It is well 
known that when the cavities are initially in photon number states a pattern of 
interference is observed at the detectors |M0lmer97al |M0lmer97b| . This is explored 



in depth in Sec. 12. II following an analytic approach first set out in |Sandersn3] . The 
approach exploits the properties of Glauber coherent states which provide a math- 
ematically convenient basis for analysing the process of localisation of the relative 
optical phase which plays a central role in the emerging interference phenomena. 
For every run of the procedure the relative phase localises rapidly with successive 
photon detections. A scalar function is identified for the relative phase distribution, 
and its asymptotic behaviour is explained. The probabilities for all possible mea- 
surement outcomes a given time after the start are computed, and it is found that 
no particular value of the localised relative phase is strongly preferred. In the case 
of an ideal apparatus, the symmetries of the setup lead to the evolution of what is 
termed here a "relational Schrodinger cat" state, which has components localised at 
two values of the relative phase. However, if there are instabilities or asymmetries 
in the system, such as a small frequency difference between the cavity modes, the 
modes always localise to a single value of the relative phase. The robustness of 
these states localised at one value under processes which obey the photon number 
superselection rule is explained. 

In Sec. 12.21 a visibility is introduced so as to provide a rigorous and operational 
definition of the degree of localisation of the relative phase. The case of initial mixed 
states is treated in Sec. 12. 3| looking specifically at the examples of Poissonian initial 
states and thermal initial states. Surprisingly it is found that the localisation can 
be as sharp as for the previous pure state example, and proceeds on the same rapid 
time scale (though the localisation for initial thermal states is slower than for initial 
Poissonian states). Differently from the pure state case where the optical modes 
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evolve to highly entangled states, in these mixed state examples the state of the 
optical modes remains separable throughout the interference procedure. Localisation 
at two values of the relative phase confuses the interpretation of the visibility, and 
a mathematically precise solution is given for the Poissonian case. 

Chapter [HI Advanced Topics on Localising Relative Optical Phase 

Chapter 01 extends the programme of Chapter |21 in a variety of directions. The 
evolution of the cavity modes under the canonical interference procedure can be 
expressed simply in terms of Kraus operators a ± b (where a and b are the anni- 
hilation operators for the two modes)^ corresponding to photodetection at each of 
the photocounters which monitor the output ports of the beam splitter. A formal 
derivation of these Kraus operators is presented in Sec. Kill The situation of initial 
states with very different intensities in each mode is addressed in Sec. 13.21 A key 
motivation for considering these asymmetric initial states is to shed light on the 
situation when a microscopic system is probed by a macroscopic apparatus. The 
discussion here takes the example of initial optical Poissonian states. The relative 
phase distributions and the probabilities for different measurement outcomes are 
found to differ substantially compared to the case of initial Poissonian states with 
equal intensities for the modes. Special attention is paid to the questions of whether 
there are preferred values for the localised relative phase and the speed of the lo- 
calisation. It is shown that the relative phase localises more slowly when the initial 
states are highly asymmetric. The transitive properties of the localisation process 
are clarified in Sec. 13. 3| again taking the example of initial Poissonian states. The 
localisation process acts largely independently of prior phase correlations with ex- 
ternal systems (although the asymmetric depletion of population with respect to 
external modes has some effect). The localised quantum relative phases have the 
same transitive properties as classical relative phases. Loss of a mode does not alter 
the phase correlations between the systems which remain. 

Sec. 13.41 explores how the canonical interference procedure could be used to en- 
gineer large photon number states using linear optics, classical feed forward, and a 

''if there is an additional fixed phase shift ^ in the apparatus, the measurement operators take 
the form a ± e*^6. The characteristic locahsation is the same. 
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source of single photons. Single photons can be "added" probabilistically, by com- 
bining them at a beam splitter and measuring one of the output ports, yielding a 
2-photon Fock state with probability 0.5 (as suggested by the well-known Hong, Ou 
and Mandel dip experiment |Hong87| ). In principle this procedure could be iterated 
to yield progressively larger number states. However, this is highly inefficient and a 
protocol is presented here which greatly improves the success probabilities. Simply 
stated, the idea is to sacrifice a small number of photons from the input states be- 
fore each "addition", in order to (partially) localise the relative phase, and then to 
adjust the phase difference to or tt. Subsequent combination at a beam splitter, 
and measurement of the output port with least intensity, yields a large Fock state at 
the other output port with much improved probability. Proposals for Heisenberg- 
limited interferometry provide one direct application of large photon number states. 
In addition, when used as the initial states for the canonical interference procedure, 
large photon number states can be used to make relational Schodinger cat states, 
and these also have potential applications. For example, when the components of 
a relational Schrodinger cat state have relative phases different by approximately 
TT, the cat state can easily be converted into a "NOON" state^ (a superposition 
of Fock states of the form |A^)|0) + |0)|A^)). NOON states are currently attracting 
considerable research interest. 

Sec. IH. 51 discussions a relational perspective on the topic of superselection rules. 
Aharonov and Susskind suggest in |Aharonov67j that superpositions forbidden in 
the conventional approach of algebraic quantum field theory can, in fact, be ob- 
served in a fully operational sense, by preparing the apparatus in certain special 
states. However, the states suggested by Aharonov and Susskind are not easy to 
prepare. Here optical mixed states with well localised relative phases and large 
intensities are presented as alternatives. These can be readily prepared and serve 
the same purpose with no loss due to the lack of purity. Finally Sec. IH.fjl suggests 
possible future extensions of the work in Chapters |21 and El for example clarifying 
the consequences for the canonical interference procedure of detector inefficiencies 

^More specifically, a phase shifter and a beam splitter can be used to convert the relational 
cat state into an "approximate NOON state" — for which the amplitude of the NOON state 
component is very much larger than other contributions. 
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tj^ical for experiments in the optical regime. 

Chapter |3l Interfering Independently Prepared Bose-Einstein 
Condensates and Localisation of the Relative Atomic Phase 

Chapter|3]looks at the interference of two, independently prepared, Bose-Einstein 
condensates which are released from their traps and imaged while falling, as they 
expand and overlap. High contrast patterns of interference have been observed 
experimentally |Andrews97j . In many theoretical treatments of Bose-Einstein con- 
densation every condensate is assigned a macroscopic wavefunction. This presumes 
an a priori symmetry breaking that endows a condensate with a definite absolute 
phase. Interference is trivially predicted on this basis. However this description 
poses various conceptual difficulties. In particular, it implies coherences between 
different atom numbers at odds with conservation of atom number, is commonly 
justiffed in terms of symmetry breaking ffelds with no clear physical relevance, and 
invokes absolute phases with values which cannot be measured even in principle. 
This is discussed more in Sec. 14.11 Such assumptions are, however, not necessary 
to predict the spatial interference of independently prepared condensates, as was 
ffrst demonstrated in detail by Javanainen and Yoo |Javanainen96bj . They studied 
the problem numerically for the case when the condensates are initially in number 
states. 

A new analysis of the interference process is presented in Sec. 14. 2| based on 
the same measurement model as used by Javanainen and Yoo. Localisation of the 
relative atomic phase plays a key role. The process of localisation is the same as 
discussed in Chapter EJfor localising relative phase between optical modes, in the case 
of an asymmetry or instability in the apparatus causing random phase shifts between 
photon detections. The visibility, defined in Sec. 12.21 of Chapter |21 to quantify the 
localisation of relative optical phase, can be translated into the current context. It 
has a simple interpretation in terms of the probability distribution for single atom 
detection. The case of initial Poissonian states is analysed in detail. The localisation 
is messy and a novel method is presented to characterise it. It is predicted that the 
relative phase distribution after the first few atom detections, defined in terms of a 
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basis of coherent states, takes the form of a Gaussian with width between 2/\/D and 
2\/2/\/D (where D denotes the total number of such detections). In particular, the 
relative phase is predicted to localise rapidly to one value, and very much faster than 
the emergence of the clearly defined patterns of interference simulated by Javanainen 
and Yoo, and others. Numerical simulations produce results sustantially consistent 
with this analysis, although the analytically derived rate of localisation is found to 
represent a slight underestimate. The discussion proceeds to open questions asking 
what, in principle, is lost when the spatial interference is analysed on the basis of 
a naive prior symmetry breaking. Finally, Sec. 14.31 reviews a recent experiment and 
several different theoretical proposals, concerning non-destructive measurements of 
the relative phases between condensates by optical means. 



Chapter HI Joint Scattering off Delocalised Particles and Localising 

Relative Positions 

Chapter El looks at localising relative positions between massive particles scatter- 
ing light. The starting point is a recent article |R,aun8j which examines two simple 
models of scattering, which are reviewed at length in Sec. 15.11 In the first "rubber 
cavity" model, a succession of photons pass through a Mach-Zehnder interferometer 
and are detected at photocounters monitoring the output ports, localising the rela- 
tive position between two delocalised mirrors in the interferometer. In the second 
"free particle" model, plane wave photons are scattered off two particles delocalised 
in a one dimensional region and are detected in the far field. The photons are either 
deflected at a definite angle or continue in the forward direction. The localisation in 
the "rubber cavity" model resembles that discussed in Chapter |21 concerning relative 
optical phase. Differently, when the light source is monochromatic the localisation 
of the relative position is periodic on the order of the wavelength of the light. For the 
"free particle" model the greatest difference is that the momentum kick imparted 
for each photodetection is variable, leading to localisation at a single value. 

The localisation in the free particle model of scattering is explored in detail in 
Sec. 15.21 and Sec. 15.31 making different assumptions. The initial states in |Rau03j 



are momentum eigenstates which are not particularly realistic. Instead the initial 
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state of the particles is taken to be thermal. The situation considered is that of an 
observer viewing a distant light source. The incident light is either forward scattered 
by the particles into the field of view of the observer or deflected, in which case the 
light source is observed to dim. Results are presented for the cases of the incident 
light being monochromatic and thermal. In both cases the localisation is only partial 
even after many detections, in contrast to the sharp localisation reported in |Rau03j . 
Possible future calculations are suggested in Sec. 15.41 

Chapter Outlook 

The Outlook suggests several possible directions for future research on this thesis 
topic. The "modus operandi" developed in Chapters El through to El can easily be 
adapted to answer a range of further questions, and can also be applied to other 
physical systems. For example, it should be possible to analyse processes localising 
the relative angle between two spin systems along the lines of calculations in this 
thesis. In another direction, it is expected that the discussion in Chapter |31 con- 
cerning the interference of atomic Bose-Einstein condensates, is relevant to systems 
of superconductors. This is suggested by the fact that Bose condensation of Cooper 
pairs — weakly bound electron pairs — plays a central role in the Bardeen-Cooper- 
Schrieffer theory of superconductivity. One concrete system where localisation of 
a relative superconducting order parameter could be investigated is that of bulk 
superconductors placed close together, and coherently coupled by a mechanically 
oscillating superconducting grain. Theoretical studies have been published which 
provide a well characterised model for the dynamical evolution of such a setup. 

Chapters El El and El which focus in different ways on the localisation of rela- 
tive quantum phases, are relevant to the debate concerning different mathematical 
characterisation of phase measurements in quantum mechanics. One possible future 
project might study the expected decorrelation of some relative number variable, 
that would be expected to accompany the localisation of a given relative phase vari- 
able. For example, what happens when the respective modes have different intensi- 
ties? This might shed some light on the dictum "number and phase are canonically 
conjugate quantum variables". In another direction, processes of localisation are 
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relevant to synchronising "quantum clocks", the subject of an extensive hterature 
concerned with problems of time in quantum mechanics. Finally, states with a well 
defined relative correlation, of the type whose preparation is discussed in this thesis, 
have potential application as pointer states in the theory of decoherence. These 
states exhibit quasi-classical properties and, in many cases are predicted to be long 
lived with respect to coupling to an environment. 
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Localising Relative Optical Phase 



This chapter looks at the interference of two, fixed frequency, optical modes which 
have no prior phase correlation, presenting and extending work published in |Cablen5] . 
A simple setup is considered wherein light leaks out of two separate cavities onto a 
beam splitter whose output ports are monitored by photodetectors. It is well known 
that when the cavities are initially in photon number states a pattern of interfer- 
ence is observed at the detectors. The dynamical localisation of the relative optical 
phase plays a key role in this process. Key studies are presented in |M0lmer97a| 
|M0lmer97bl ISanders03] and some results are also reported in |Chough97| , all fo- 
cusing on the case of initial number states with the same photon number in each 
mode. Some results for the first and second photodetections are also presented in 
|Pegg05| for the general case of mixed initial states with zero optical coherences. The 
work here adopts the analytic approach introduced in |Sandersn3] and represents a 
substantial development of the topic. 

Sec. 12.11 looks in detail at the case when the cavity modes are initially in pho- 
ton number states. It first explains the basic setup and provides some intuition as 
to why interference phenomenon are observed. M0lmer's numerical treatment of 
the problem in |M0lmer97at |M0lmer97b| is also summarised. Adopting a basis of 
Glauber coherent states. Sec. 12.1.11 investigates the evolution of the relative phase 
distribution for all possible measurement sequences. Expressions for the asymptotic 
forms of the relative phase distributions are also given. The probabilities for all pos- 
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sible measurement outcomes a given time after the start are evaluated in Sec. I2.1.2| 

and it is found that no particular value of the localised relative phase is preferred 

in this example.^ In addition it is seen that an additional fixed phase shift in one 

of the arms of the apparatus does not alter the experiment. The emergence in 

an ideal experiment (one without phase instabilities) of what is termed a relational 

Schrodinger cat state is discussed in Sec. l2.1.Jl Sec. 12 . 1 .41 discusses the robustness of 

states sharply localised at one value of the relative phase, including a brief summary 

of the relevant simulations in M0lmer97b| . 

A visibility is introduced in Sec. 12. 21 as a means to rigorously quantify the degree 

of localisation of the relative phase. It ranges from (no phase correlation) to 1 

(perfect phase correlation). In Sec. 12.31 the analysis of Sec. 12.11 is extended to the 

case of mixed initial states, which is more realistic experimentally. Specifically the 

case of Poissonian initial states is treated in Sec. 12.3.11 and of thermal initial states 

in Sec. 12.3.21 Differently from the pure state case where the localised state of the 

two cavity modes is highly entangled, in these mixed state examples the state of 

the cavity modes remains separable throughout the interference procedure. The 

visibilities for the final states are computed for all possible measurement sequences. 

Surprisingly it is found for both examples that the localisation can be as sharp as for 

initial pure states, and proceeds on the same rapid time scale. In fact the localisation 

turns out to be slightly slower for the thermal case. For Poissonian initial states the 

visibility jumps from to 1/2 after just one measurement, whereas for thermal initial 

states is jumps to 1/3. When the interference procedure involves the detection of 

more than one photons, and photons are registered at both detectors, the relative 

phase localises at two values. The visibilities for these cases underestimate of the true 

degree of localisation, and a solution is presented for the Poissonian case. Finally 

12. 3. 31 discusses the consequences of errors from the photodetectors, which is a feature 

of any real experiment. 

^Different outcomes at the detectors are associated with locahsation at different values of the 
relative phase. Some measurement outcomes are more likely than others. However the least 
likely measurement outcomes cause localisation at values which are closer together. Overall the 
localisation process does not "prefer" any particular range for the localised relative phases for this 
choice of initial states. 
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2.1 Analysis of the canonical interference proce- 
dure for pure initial states 

A simple operational procedure for both causing and probing the localisation be- 
tween two optical modes is depicted in Fig. 12. 1[ Two cavities containing and 
M photons respectively at the start (and thus described by pure initial states |A^) 
and |M)) both leak out one end mirror (via linear mode coupling). Their outputs 
are combined on a 50 : 50 beam splitter, after which they are detected at two 
photocounters. 



/^Cavity A 



Cavity B 7\ 

mA 

Phase 
shifter 




Photodetector Photodetector Z)^ 

Figure 2.1: Photon number states leak out of their cavities and are combined on a 
50:50 beam splitter. The two output ports are monitored by photodetectors. In the first 
instance the variable phase shift ^ is fixed at for the duration of the procedure. 



Despite the cavities initially being in Fock states with no well-defined relative 
phase it is well known that an interference pattern is observed at the two detectors. 
The interference pattern can be observed in time if the two cavities are populated by 
photons of slightly differing frequencies or, as in standard interferometry, by varying 
a phase shifter placed in one of the beam splitter ports. Despite the evolution for 
the system taking place under an effective superselection rule for photon number, 
coherence phenomena depending on the conjugate phases are thus observed. The 
reason for this contradiction with the dictum "number and phase are conjugate 
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quantities" may be understood as follows. 

Consider the case after a single photon has been detected at one of the detectors. 
Then the new state of the two cavities is ^jj^\N-l)\M) ± ^jj^\N)\M-l) i.e. 
it is entangled. It is simple to show that the second photon is much more likely to be 
registered at the same detector. The exact ratio of the probabilities of being counted 
at the same detector and at the other is N'^+M'^-N-M+ANM to N'^+M'^-N-M. 
When N = M this ratio is strictly greater than 3, and tends sharply to infinity as 

and M approach 1. This is in agreement with the phenomenon demonstrated 
by the well-known Hong, Ou and Mandel dip experiment |Hong87| , whereby two 
uncorrelated and identical photons, simultaneously incident on the input ports of 
a 50 : 50 beam splitter must both be registered at the same output port. Further 
detections lead to a more and more entangled state. It is not so surprising then 
that detections on an entangled state lead to some form of interference pattern. In 
essence after a small number of detections the relative number of photons in each 
cavity is no longer well defined, and so a well defined relative phase can emerge. 
Note that this is only possible if the beam splitter, detectors and cavities all have 
well-defined relative positions. 

One method for confirming this intuition is to use a quantum jumps approach 
(for a review of quantum jump methods see |Plenio98] and references therein), and 
numerically simulate such a system through a number of detection procedures. Such 
an approach was taken in Sec. Ill of |M0lmer97a| and Sec. 4 of |M0lmer97b| , and 
some important features of those studies are highlighted here. M0lmer assumes a 
frequency difference ujb—uja between the two cavities, which gives rise to an interfer- 
ence pattern in time. The cavities are assigned equal decay rates F for leakage onto 
the beam splitter, leading to exponential decay laws for the cavity intensities and the 
photon count rates. M0lmer looks at the case that both cavities start with the same 
number of photons, say, and the chosen values of are large. In |M0lmer97a| the 
choice of parameters is A^ = 10^ and Uh — Ua = 1000F, and in |M0lmer97b| N = 5000 



and cufo — = 30F.^ Sharp localisation of the relative phase occurs as the first few 
photons are detected, which manifests itself as continuous oscillation in the count 



^In contrast the analysis in this chapter for pure initial states applies for arbitrary initial cavity 
photon number N. 
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rates at the two photodetectors as further photons are registered. This oscillation 
takes place on a time scale l/{uJh — uja), with a high count rate alternating between 
the detectors. M0lmer also looks at the evolution of the quantity Q = {ilj\a^b\i/j) 
(where is the state of the two cavities and, and b are creation and annihilation 
operators for cavities A and B respectively).'^ This quantity determines the count 
rates at the detectors and makes a rapid transition from a random to a harmonic 
evolution as the first few photons are registered. However simulations such as these 
yield little in the way of physical insight. The following discussion is based instead 
on an analysis introduced in |Sanders03] . 

2.1.1 Evolution of the localising scalar function 

To begin, the initial state lipi) = |A^)|M) of the cavities is expanded in terms of 
coherent states \a), 

\^j)=^ rY'rf^d0e-*(^^+*-^*)|a)|/?) (2.1) 
Jo Jo 



with a = ^/Ne'^, f3 = y/Me''*', and the normalisation Af = l/A/nAr(A^)nM(M)47r2 
where n„(/i) = /i"e~'^/n! is the Poissonian distribution. For the properties of coher- 
ent states refer |Glauber63j (for pedagogical reviews see |Klauder85l Gerry04| ) . At 



first the normalisation in Eq. ()2.H1 will be ignored for simplicity. 

Consider now the case that a single photon is detected at either the left detector 
Dl or the right one Dj^. Since only the change of the state of the cavity modes is 
of interest the exterior modes may be treated as ancillae, and corresponding Kraus 
operators Kl and Kn describing the effect of the detection on the cavity modes only 
may be determined (for an explanation of the basic properties of Kraus operators 
see for example Chapter 8 of [ NielsenOO] ) . Treating in a basis of coherent states the 
leakage of the cavity populations into the ancillae, and subsequent combination at 
a beam splitter and photodetection, as is done in jSandersOS] , immediately suggests 
the form of Kl and K^i as proportional to a±6 (where a and b are annihilation oper- 
ators for the modes in cavity A and B respectively). The constant of proportionality 

•^More specifically he considers 2Ke{Q) / {2N — D) (where Re(-) denotes the real part of some 
complex number and D denotes the total number of detected photons). 
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depends on the transmittivity of the cavity end mirrors. For the purposes here of 
examining the main features of the locahsing relative phase, this resuh is assumed 
to be correct, and a careful verification is delayed to Sec. l3.1l of Chapter El'' It may 
be observed immediately that these simple expressions for Kl and constitute a 
complete measurement process. This reflects the assumptions that the beam splitter 
is lossless and acts as a unitary process on the ancillae. The equal weighting of a 
and b reflects the assumptions that the beam splitter is 50 : 50 splitting, and that 
rates of the leakage from both cavities onto the beam splitter are the same. 

In the event that some number I of photons are registered ai and r at Dji, 
the state of the two cavities evolves as follows: 

oc K[K'j,jjded<l)e~''^'''+''*^\a)\(3) 

oc j l'ded(f)e-'^^'+^'^\a - (3y{a + (3y\a)\(3) 

The scalar function, 

Q^r^ia-pYia + py, (2.2) 

encodes information about the localisation in the relative phase which occurs be- 
tween the two cavities. It should be pointed out that in this example there is some 
ambiguity over the definition of Ci^r- In expanding a Fock state in terms of coherent 
states, |A^) = ^ ^ J^^ ^e~'^'^'^\^/me^'^) , the number m is free to take any positive 
value — the key feature of the expansion for a state with photons is that the 
integral must encircle the origin in phase space times (where the phase space is 
the complex plane for the variable labelling the basis of coherent states). Altering 
the relative amplitudes of a and (3 in Eq. ()2.1|) will alter the relative phase distri- 
bution contained in Ci^r- When the cavities begin in the same photon number state 
the natural choice is to set the amplitudes of the coherent states the same for both 
modes. 

■^The operators a±b are also derived in |M0lmer9'7al |M0lmer97b| and |Pegg05| using different 
arguments. Detailed calculations later in this chapter do not assume this result but treat the 
evolution of the cavity and ancilla modes in full. 
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Figure 2.2: The evolution of C;^j.(0, (/>). In (a) localisation about Aq = tt after 1, 5 
and 15 counts when photons are recorded in the left photodetector only, (b) localisation 
about Aq = ib2arccos (l/-v/3) ~ 1-9 after 3, 6 and 15 counts when twice as many photons 
are recorded in the left detector as the right one. The symmetry properties of the Kraus 
operators and Kji cause Ci^r to have multiple peaks (either one or two for A ranging 
on an interval of 27r). 



For the purposes here it is sufficient to focus on the symmetric case N = M 
when both cavities begin in the same state. It should be noted however that the 
physics of the highly asjTiimetric case is somewhat different (for further discussion 
refer Sec. IH.2l of ChapterEI). Setting |a| = \(3\, the "localising scalar function" Ci^r 
takes the form, 



i<f>\r 



(2.3) 



(4iV) 



a+O/2(_,);gia-f.)(e+0)/2sin«Acos^f 



where A = (j) — 6. Factors that do not depend on 6, cf) will be ignored for the moment, 
since they will be taken care of by normalisation. 

Of particular interest is the behaviour of C;,r(^; 0) as the total number of detec- 
tions / + r gets larger. Asymptotic expansions |R,owenij for Ci^r{0, 0) can be used to 
examine this limit. When photons are detected at both detectors. 



, / A r A I / t r 

sm' fcos*^ f ~ \; , M 

2 2 1 y /_^^i-fr 



I' 



exp 



^(A-Ao)' 



(2.4) 



where Aq = 2 arccos ^/r/{r + 1) when A takes values between and tt, and Aq = 
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27r — 2 arccos ^/r/{r + I) between vr and 27r. Aq denotes the values of the relative 
phase around which the localisation occurs. When all the photons are detected at 
one detector the appropriate expressions are 



cos 



sm 



r A I 
2 I 

I A I 



exp 
exp 



-TT, TT 



(A- 



-vr 



for A G ^ 

for A G [0,27r]. 



(2.5) 



In every case asymptotically \Ci^r\ takes the form of a Gaussian distribution with 
width (two standard deviations), A/y/r (or 4/v^) when the counts are all at one 
detector or 2\^/ y/l + r otherwise, decreasing with the total number of detections. 

As / + r grows the state of the two cavities evolves into a superposition (over 
global phase) of coherent states with an increasingly sharply defined relative phase. 
A plot showing the evolution of Ci^r{(^,4') is shown in Fig. 12.21 This localisation 
in the relative phase is responsible for the interference phenomena seen at the two 
detectors, as was examined numerically in |M0lmer97at |M0lmer97b . 



2.1.2 The probabilities for different measurement outcomes 

The value Aq at which the relative phase localisation occurs depends on (the ratio 
of) the specific number of photons / and r detected at each detector. This is of 
course probabilistic. Letting Pi^r denote the probability of detecting I and r photons 
at the left and right detectors respectively, a complete expression for Pi r may be 
obtained by a simple heuristic treatment of the djTiamics as in |Sanders03] . It is 
supposed that population leaks out of each cavity modes into an ancilla according 
to a linear coupling with parameter e, where e is small. The ancillae evolve under 
the action of a 50 : 50 beam splitter. Subsequent action of the projection operators 
|/)(/| and \r){r\ on the ancillae corresponds to the detection of / photons at the left 
detector and r photons at the right detector. 

It is helpful here to take a brief diversion to clarify the mathematics of a general 
dynamical evolution under some linear mode coupling, which is supposed here to 
govern the leakage from each cavity, and the action of the beam splitter. Denoting 
the creation and annihilation operators for some pair of modes c^, c and ,d, a linear 
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mode coupling Hamiltonian is a linear combination of energy conserving terms of 
the form, 

where and ^ are real parameters. Evolution over some time t is then given by the 
unitary operator, 

C^imc iv, = exp (^-iHimct^ 



= exp I (pe 



where (/? = M. One particularly important property of linear mode couplings is that 
they evolve products of coherent states to products of coherent states. Specifically 
for coherent states with complex parameters a and 

Uimc (^^,0 l«) \P) = |Q;cos^? + /3e"*^sin^?) \-ae'^ siwd + (3 cos-d) . 

Leakage into the vacuum is described hy P — and a small value for i? so that 
cos'j? = \^l — € and sini? = ^/€ where e is small (and say ^ — n). The hnear mode 
coupling then acts as. 



|a)|0) — > \Vl^a)\^/ea). 



For a 50 : 50 beam splitter i? = 7r/4 so that (with the additional phase shift ^ = 
say), 

\a) \3) ■ 



a + (3 



V2 



V2 

Many calculations are facilitated by considering how a linear mode coupling trans- 
forms the creation and annihiliation operators. This is given by,^ 

U\e,(t))cU {9,4)) = cos + 6-^"^ sin 
U^{e,(j))dU {e,(j)) = -e'^smOc + cosed. 



"''Sometimes c and d arc referred as "input" field operators and denoted Cin and din say, while 
after unitary transformation the operators are called "output" field operators so that Cout = U^Cinll 
and dout = U^dinU. This notation is not used here. 
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Returning to the interference procedure the full expression for the cavity modes 
after the measurement process has acted is, 

AfJd9d(f)e-'^^+''^Ci,r{e,e,^) \^/T^a) \VT^P) , 

where the normahsation factor J\f — (njv(-/V)47r^)~^ and, 



-a + (3 



-a + (3 



C,.(6, 0, 4>) = {r Ve^j [I ^'^yf^/ ' (2-6) 
where (r| and (/| denote photon number states. The probability Pi^r is given by. 



X Ci,r{e,e,(l))Ci,r{e,e', <!)')* 
X {VT^a'\\VT^a){y/T^(3'\\VT^(3). 

To simplify this expression one can at first proceed naively, taking the overlaps of 
the basis states to be zero, to derive an approximate expression for Pi^r which can be 
compared numerically to the real values. It might be guessed that the approximate 
expression thus derived will hold good when several, but not too many, photons 
have been recorded so that C; is narrow while the amplitudes ^/l — ea, \/l — e/3 
are still large. Hence assuming, 

{a'\a) = exp (— |a — a'p) ~ 5(0 — 0') 
W) = exp(-|/3-/?f)~<5(^-^'), 

and using the relation for the gamma function r(.), 

^"'•2'' ^ r(r + 0.5)r(/ + 0.5) 



-^|^cos^'^#sin2'- 



lo Jo ' ' 7rr(r + / + l) 

the following approximation for Pi^r is obtained. 



l,r ' 



{r + l)\ 



;r + /)!r(r + 0.5)r(/ + 0.5) 



rill 7rr(r + / + l) 



(2.7) 
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Figure 2.3: A plot of the exact values of the probabilities Pi^r for all the possible mea- 
surement outcomes to the procedure a finite time after the start, against the absolute 
value of the relative phase which is evolved. The initial state is |20) |20) and the leakage 
parameter e, corresponding roughly to the time, has a value of 0.2. Each spot corre- 
sponds to a different measurement outcome with / and r counts at detectors Di and ir- 
respectively. The value Aq of the relative phase which evolves in each case is given by 
2 arccos ( Y^r/(r + /) ) . 



It is found numerically that this approximation is surprisingly good, and applies 
quite generally whenever e is small. The fractional error goes roughly as 0.6e, growing 
linearly with the leakage parameter. As for its general features Pi^r is seen to be a 
product of a global Poissonian distribution in the total number of detected photons 
/ + r and a second function depending on the precise ratio of counts at Di and Dr. 

A plot of the exact values for the probabilities for different measurement 
outcomes is plotted in Fig. 12.31 for typical parameter values e = 0.2 and initial state 
1 20) 1 20), with each spot corresponding to a possible outcome, e corresponds to a 
time parameter, an approximation which holds good provided e is not too large. 
The distribution Pi ,, gives the likely degree of localisation of the relative phase A 
a finite time after the start of the procedure, and the values Aq which are picked 
out. Looking at the precise distribution in Fig. 12. 3| it is seen that given e = 0.2 it is 
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most likely that 7 photons (approximately 2eN) have been counted (corresponding 
to the ridge). For the most probable outcomes all the photons are counted at one 
detector, in which case the relative phase localises at or tt. However the density 
of points is greatest about Aq = |, and for these outcomes there are approximately 
equal counts at both detectors. Overall no particular value of the localised relative 
phase is preferred in this example. 

Finally it is interesting to ask how the probabilities for different sequences of 
detections are altered if the phase shift ^ between cavity B and the beam splitter 
(refer Fig. 12.11) is a fixed at a value not equal to for the duration of the procedure. 
The phase shift alters the function Ci^r{^, 0)? Eq. ()2.6|) . according to a translation 
in the relative phase i.e. A = — 0— i>A + ^. The Kraus operators Ki and Kj. for 
the altered apparatus are a ± e*^6. Looking now at the effect of r "right" detections 
and 1 "left" detections on the initial state |A^)|A^) in a basis of Fock states, 

\N)\N) {a + e%y{a-e%)^\N)\N) 

= 4t'orf,e^('+'-^-)«|iV-j)|iV-/-r+j), (2.8) 

where the dj are some coefficients independent of ^, reveals that the dependence on 
C, drops out on taking the norm of the final state. Hence the probabilities P; are not 
affected by the additional phase shift despite the asymmetry in the apparatus. The 
phase shift ^ is seen to be relevant operationally only if altered at different stages 
during the procedure. 

2.1.3 Symmetries of the localising procedure and relational 
Schrodinger cat states 

Once a given measurement outcome has occurred with / and r counts at the left and 
right detectors respectively the resultant state of the two cavities has two symme- 
tries. This is demonstrated by the explicit form of Ci^r, Eq. ()2.3|1 . and in Fig. 12.21 A 
271 translational symmetry identifies physically identical phases. In addition there 
is symmetry in C/ r about A = 0. This exists because the procedure as described so 
far localises the absolute value of the relative phase. 
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When photons are detected at both ports Ci^r is peaked at two different values 
±Ao. Looking at the asymptotic form of Ci^r as I and r tend to large values it is 
observed that the state that emerges, \ipoo), takes the following form: 



g-i|7|2Ao|^giAo^^gi|7pAo|^g-iAoJ ^ (2.9) 



where I7) = | |7|e*^) and I7I = y — (^) • The relative component of the two mode 
state, contained in the square brackets, is a superposition of two coherent states with 
the same amplitude but different phases ±Ao — ordinarily called a Schrodinger cat 
state, l^/'oo) has in addition a sum over all values of the global phase 9. A state of the 
form Eq. ()2.9p could be termed a relational Schrodinger cat state. It may be asked 
why values of the relative phase with the same magnitude but opposite sign may not 
be identified as equivalent, in particular since either one of the two components of 
Eq. ()2.9|) separately gives rise to the same probabilities for further detections at the 
left and right photocounters. As an example of why the two values of the relative 
phase are not equivalent consider the case that a relational Schrodinger cat state is 
formed with |Ao| — 7r/2, and the phase shifter (refer back to Fig. \2.\\ initially fixed 
at ^ = 0) is subsequently adjusted by 7r/2. With high probability later photons will 
all be detected at one detector, which is randomly the "left" or "right" one varying 
from run to run of the experiment. This behaviour is not consistent with a state 
sharply localised at one value of the relative phase. 

Creating the superposition Eq. ()2.9|) would however be experimentally challeng- 
ing as it requires perfect phase stability. In practise it is found that the relational 
Schrodinger cat is sensitive to any asymmetry or instability in the system. The effect 
of a randomly varying phase is to cause localisation about one particular value of 
the relative phase. This phenomenon is evident in the numerical studies of M0lmer 
|M0lmer97at |M0lmer97b| . These incorporate a slight frequency difference between 
the two cavity modes causing the free evolution to have an additional detuning term 
expi(co'fe — uja)i>^ht. Combined with the random intervals between detections, this 
means that the process can be described by Kraus operators a ± e*'^6 where the 
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phase a takes random values for each photodetection. The relative phase then takes 
a unique value varying randomly for each run. A dynamically equivalent process 
occurs when atoms from two overlapping Bose Einstein condensates drop onto an 
array of detectors and are detected at random positions; a detailed discussion of this 
point is given in Sec. 14.2. fl of Chapter HI 

In the case of an idealised setup for which the phase shifts throughout the ap- 
paratus remain fixed, one component of the relational Schrodinger cat state can 
be removed manually. Suppose that after / and r photons have been detected at 
Di and Dr in the usual way, the phase shifter is adjusted by ±Ao, and then the 
experiment is continued until a small number of additional photons have been de- 
tected. The phase shift translates the interference pattern in such a way that the 
additional counts will occur at one detector (with high probability). The additional 
measurements eliminate the unwanted component of the cat state and confirm a 
well-defined relative phase. 

2.1.4 Robustness of the localised states and operational equiv- 
alence to tensor products of coherent states 

The next important feature of the canonical interference process concerns the ro- 
bustness of the localisation. In the limit of a large number of detections, the state 
of the two cavities becomes equivalent to 



with I7) = ||7|e*^) some coherent state. The coherent states, being minimum uncer- 
tainty gaussian states, are the most classical of any quantum states (see for example 



bust. However lipoo) is a superposition over such states, and this could potentially 

affect the robustness. That this is not the case can be understood by noting that the 

superposition in Eq. ()2.10|1 is summed over the global phase^ 6 of the coherent states. 

® "global phase" here is not referring to the always insignificant total phase of a wavefunction, 
but rather the phase generated by translations in photon number: e"^". This is still a relative 
phase between different states in the Fock state expansion of a coherent state. 




(2.10) 



Chapter 2 of |Scully97l ). Thus states of the form |7)|7e*^°) are expected to be ro- 
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Under evolutions obeying an additive conservation of energy rule (photon-number 
superselection), which is essentially the extremely good rotating- wave approxima- 
tion of quantum optics, this global phase becomes operationally insignificant (refer 
|Sandersn3] ). This is discussed further in Sec. 13 .^^ Chapter IHl 

Evidence of the robustness of the localised state Eq. ()2.1()j) has also been pro- 
vided by numerical studies reported in Sec. 5 of |M0lmer97b| , and a summary of 



M0lmer's findings is as follows. The setup considered was similar to that depicted 
in Fig. 12.11 except the cavities had a slight frequency difference, and each cavity 
was assumed to be coupled to an additional, independent, reservoir. M0lmer first 
looked at the scenario in which the reservoirs acted as additional decay channels 
and had decay rates equal to that which determined the leakage onto the beam 
splitter. The characteristic oscillation of the counts rates at the two measured ports 
of the beam splitter was found to be unchanged by this additional leakage, although 
the intensities of the cavity fields and the counts rates decayed twice as fast. In 
the second scenario considered the reservoirs were thermal and photons could feed 
into the cavity as well as out. The photodetector count rates were more difficult 
to interpret in this example. The diagnostic variable 2'Re{{ilj\a^b\tp)), where a and 
b denote annihilation operators for the two cavity modes and the state of the 
cavities, was observed to oscillate consistent with localisation of the relative phase 
within the first few detections. However the oscillations were not as smooth as in the 
absense of the reservoirs. Furthermore the oscillations were disrupted and partially 
reestablished, a phenomenon which was found to coincide with the cavities being 
nearly empty. It was concluded that the relative phase correlation of the field modes 
is not robust when the cavities are populated by only a few photons. 

Finally a state of the form Eq. ()2.10j) is, for any processes involving relative 
phases between the cavities, operationally equivalent to a tensor product of pure 
coherent states for each cavity |7)|7e*^''). However, because of the phase factor 
g-2«|7|26»^ the state is in fact highly entangled. Expanding in the (orthogonal) Fock 
bases (as opposed to the non-orthogonal coherent states) the state is seen to be of 
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the form: 



1^00) = j fe-^'l^l^^|7)|7e^^°> 
00 „ 



n,m=0 
2|7P 



Evn2|.|^^(l7l^n„(|7|2) e^-^°|2|7|'-m,m> (2.11) 



m=0 

where |n,m) denotes a product of photon number states, n.(.) denotes a Poissonian 
factor and 2|7p is a whole number of photons. Note that every term in the su- 
perposition has the same total photon number, consistent with a photon number 
superselection rule. 

2.2 Quantifying the degree of localisation of the 
relative phase 

A definition of the visibility suitable for rigorously quantifying the degree of localisa- 
tion of the relative phase for a prepared two mode state is illustrated by Fig. 12.41 It 
is supposed that the second mode undergoes a phase shift r before being completely 
combined with the first at a 50:50 beam splitter. The expected photon number at 
the left port is then denoted I{t). This intensity is evaluated for all possible phase 
shifts r, allowing a visibility for the two mode optical state to be defined in terms 
of the difference between the maximum and minimum values as follows, 

^ (-^max -^min) / (-^max ~l~ -^min)' (^•-^^) 

By definition the visibility takes values between and 1. 

For a product of photon number states |A^)|M) the action of a phase shifter 
on the second mode merely introduces an irrelevant factor of e*^^, and hence the 
intensity /(r) is constant for different phase shifts r, and the visibility \^ is 0. 
In a similar way the visibility is for any product of mixed states diagonal in 
the photon number basis, such as the product of Poissonian states in Eq. ()2.13j] 
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Mode 1 



Mode 2 



Phase 
shifter 




Photodetector 

Figure 2.4: /(r) is the intensity at the left output port after the second mode undergoes 
a phase shift of r and is combined with the first at a 50 : 50 beam splitter. This intensity 
is evaluated for all possible settings of the phase shifter. Extremising over r, the visibility 
for the two mode state is defined as V = (/max — /min)/(/max + Imin)- 



(discussed later in Sec. 12.3. On the other hand, for a product of coherent states 
|yiVe^^)|yiVe^(^+^«)) and, |V^oo) Eq. fU^i and poo Eq. fT^ (discussed later in 
Sec. 12.3. ip summed over the global phase, and all three with exactly one value 
Aq for the localised relative phase, it is easily shown that J(r) is proportional to 
cos^ ( ^°2^^ ) • -^(t) is then maximized if the phase shifter is set to r = — Aq and for 
r = — Aq + TT. Therefore the visibility is 1 for these three examples for which the 
relative phase is perfectly correlated. 



2.3 Mixed initial states 
2.3.1 Poissonian initial states 

The example of Sec. 12. while usefully illustrating many features of relative locali- 
sation, is not experimentally accessible due to the assumption of the availability of 
large photon number, initially pure, Fock states populating the cavities. In partic- 
ular, if looking for a mechanism by which relative localisation occurs naturally in 
our interactions with surrounding objects, the previous example is somewhat im- 
plausible as it stands, in as much as it would suggest that macroscopic levels of 
entanglement are necessary to localise relative degrees of freedom. 
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With this in mind a more reahstic scenario is considered in this section. While it 
is implausible that the cavities are populated by large Fock states, it is not implau- 
sible that they are populated by a large number of photons, and that only the mean 
number N of photons is known. In such a situation the quantum state of the cavity 
would be assigned a Poissonian distribution over photon number under a maximum 
entropy principle (refer |Jaynes03| ). Alternatively, the cavities may be populated by 
light from independent lasers, for which standard laser theory leads to the photon 
number distribution being Poissonian (refer Chapter 11 of |Scully97| ). 

The canonical localisation procedure discussed in Sec. 12.11 is now reconsidered 
assuming that the initial state of the cavities is. 



Pi = J]n„(iv)|n)(n|® J]n^(iv)|m)( 

n 
1 



47]-2 



ded(j)\a){a\®\(3){f3\, (2.13) 



where a = vN expiO and /? = vN expi(f). The evolution of pj is as follows, given 
that I and r photons are detected at the left and right detectors respectively: 

PI KiK^piKlK^; 

oc K[K]^ j jdOdcj) \a){a\ ® \(3){(3\K^Ik^^ 

oc j jdOdcj) \a + /^I^H" - /5rV)("l ® 
oc j jdOdcP |G,,(^,0)|»(a| ® |/?)(/?|, 

with Ci^r{0,(j)) as in Eq. ()2.3j) . Clearly the discussion about the localising nature of 
Ci^r{9, 0) applies equally well in this case. The expression Eq. ()2.7|) approximating 
the probabilities for different measurement records when the initial state is a product 
of Fock states is exact for a product of Poissonian states. In the limit of a large 
number of detections, 

/•27r r/n 

Poo= / — |ae*^'')(ae^^°|. (2.14) 
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Quite remarkably the relative phase localisation of the mixed states is just as sharp 
and just as rapid as that of the pure statesl 

However, in striking contrast to the large entanglement observed for the pure 
state case considered in Sec. 12.11 the state of the two cavities in this example remains 
manifestly separable (unentangled) throughout. This fact might seem surprising at 
first sight. If the calculation had been done in terms of photon number states, 
products of photon number states |A^) |M) would evolve to highly entangled 
states. However the photon number variables are summed according to Poissonian 
probability distributions and this changes things. As an example of the fact that 
a mixture of entangled states can be separable consider the case of an incoherent 
mixture of all 4 maximally entangled (Bell) states of 2 qubits (refer jNielsenOT)] ) . This 
is equivalent to the completely mixed state 1/4 which is certainly separable. Any 
two-party mixed state of the form pa ® Pb is certainly unentangled. Furthermore, 
any "convex sum" T,a,i,Ca,bPa ® Pb, where all the pa and pb are bona fide density 
operators, has only classical correlation and by the standard definition is separable. 
This is exactly the case for the mixed state calculation considered in this section 
both before and after the measurement process has acted. A calculation in terms 
of mixtures of photon number states would likely be much messier, but ultimately 
must be consistent with this. It should be noted further here that poo in Eq- ()2.14j] 
has the interesting feature of being formally separable but not locally preparable 
under a superselection rule (or equivalently lack of a suitable reference frame), a 
feature first noted in |Rudolph'oT| . 

After involved calculation, extending methods and results developed earlier in 
this chapter, a simple expression for the intensity and the visibility can be found 
for the case of two initial Poissonian states and an idealised experiment for which 
the phase shifts throughout the apparatus remain fixed (for the derivation refer 
appendix lA. II and specifically Sec. IA.1.2|l : 



/(r) oc r cos^ | + / sin^ | + 



1 

2' 



(2.15) 



r — / 



(2.16) 



r + l + 1' 
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If the detections are all at one detector, the right one say, Eq. ()2.1(jj) simplifies to 
r/ (r + 1) which tends rapidly 1, and is 1/2 after just one detection. However, it is 
also seen that the expression decreases to as the proportion of counts at the two 
detectors becomes equal. This does not refiect less localisation in these cases but 
is an artefact of the definition of the visibility. It is easy to see that if the state 
of the two cavities is localised at two values of the relative phase these will both 
contribute to the intensity at one port in the definition of the visibility; changing 
the phase shift r will tend to reduce the contribution of one while increasing that 
of the other so that overall the variation in the intensity is reduced. In more detail, 
the visibility is underestimated when r ^ / because of the occurrence of two peaks 
in Ci^r approximately vr apart. However, in realistic situations such localisation at 
multiple values is expected to be killed by small instabilities in the apparatus (as in 
the case of relational Schrodinger cat states in pure state case) and thus the visibility 
is expected to tend to 1 in all cases. 

The definition of the visibility can in fact be modified so as to provide a math- 
ematically rigorous measure of the degree of localisation when the relative phase 
parameter is peaked at more than one value. The localising scalar Ci^r{di 4') for the 
case of initial Poissonian states (with the same average photon number) is the same 
as for the case of initial Fock states (with the same number), and is illustrated in 
Fig. 12.21 In particular it should be observed that C^r{G,4') = when (p = 6. It is 
clear then that the final state of the two cavity modes may be considered a sum 
of two separate components, one localised at + Aq = 2 arccos \/r/{r + I) with the 
relative phase parameter ^ = varying on to 7r/2, and another at — Aq with ^ 
restricted to —7i/2 to 0. Hence an improved measure of the degree of localisation is 
achieved by computing the visibility for one of these components taken separately. 
This is done in Sec. IA.1.31 of the appendix lA.ll where full analytic expressions are 
provided. Given a total number D of detections at both detectors, the visibility V 
is found to be smallest when all the photons are detected at one photocounter, with 
V = D/{D + 1). The same expression is obtained from Eq. ()2.16|) when r = or 
Z = as would be expected as the localisation process picks out a unique value for 
the relative phase in these two cases. V is found to be greater when there are counts 
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at both photocounters, and in the "fastest case" where there are equal counts at 
both photocounters, 

^ ^ 2[T{D/2 + l)f 



{D+l)[r{D/2 + l/2)\ 
A- ^ 



where F (•) denotes the gamma function, A is 1.27 when D — 2, 1.13 when D — A 
and 1.05 when D — 10, decreasing rapidly to 1. 



2.3.2 Thermal initial states 

In the examples presented above there was no limit to how sharp the localisation 
could be. The case that both cavities are initially populated by thermal states 
with the same mean photon numbers N is now examined along the same lines. A 
thermal probability distribution in the photon number n with mean N is given by 
7V"/(1 + TV)""*"^, decreasing monotonically with n, and as such the initial state of 
the two cavities is,^ 

Pi — y -, =T — 7\n){n\ ® > m)(m 

= f fdndmded(l)e-^\"\'+\^\'^/''\a){a\^\(3){(3\, 

A-K^N^ J J 

where a — y/^ex-^iO and (5 — \/Wiex^i(t). Under the measurement of I and r photons 
at the left and right detectors respectively, 

PI KiK^^piK^lK^^ 

oc fd^ad^(3e<\''\"^^\")'^KiK'R\a){a\(^W 



oc 



oc 



y"dW2/3e-(l'^l'-^^l')/^|a+/3|2n«-/3ri«)(a|(8)|/3)(/3| 
/.^„.^;,e-(l"P*)/^|a„(..«.«.^)|>X„|«IW|. 



^Note that in the analysis here the thermal light is assumed to be single mode. A valuable exten- 
sion would be to consider the interference of multimode thermal light, characteristic of unfiltered 
chaotic light from natural sources. 
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As in the previous example of initial Poissonian states the state of the cavities 
remains separable throughout the localising process. 

Unlike previous examples \Ci^r {n, rh, 9,(j)) p does not provide a simple picture of 
the localisation of the relative phase due to the additional dependence on the mean 
photon number variables. One possibility would be to average the scalar function 
over the two photon number variables leading to a function, 



dndm ^-n/N^-M/N^c^^^^-^^^Q^^^^2_ (2.17) 




^0 



N N 



After normalising this can be shown to be, 

1 pco poo 

|Q^,(A)|^ = / / dndThe~''~'^\Vfi + Vrhe'^\"\-Vfi + Vrhe'^\ , 



^0 



where A = (f) — 6. \Ci^r is independent of N. However it is not necessary to 
consider such a function in detail to quantify the rate and extent of the localisation 
of the relative phase parameter. 

An intensity and a visibility can be computed as in the Poissonian case above 
(see appendix IA.2jl . For an arbitrary measurement record the results are. 



/(r) oc /cos^ f +rsin2 f + 1 (2.18) 



Vir= , ■ (2.19) 

If all the measurements occur in one detector, the right one say, the visibility is 
r/ (r + 2) which is 1/3 after just one detection and which tends to 1 rapidly - but 
slower than in the Poissonian case. Again it is found that the localisation of the 
relative phase parameter can be as sharp for initial mixed states as for initial pure 
states. In addition, the localisation for initial thermal states is rapid, but slower than 
for initial photon number states or Poissonian states. 

Averaging over all detection outcomes, an expected visibility a finite time after 
the start of the procedure T,i^rPi,rVi^r can be computed. An exact expression for the 
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Figure 2.5: Expected visibilites for (a) an initial product of two Poissonian states 
(plusses) and (b) an initial product of two thermal states (crosses), with average photon 
number N for both cavities. 



probabilities of different measurement sequences derived in (see appendix IA.2jl is, 



Pjr 



(f + eN) 



r-r\r+l+2 ■ 



(2.20) 



This notably has the form of the probabilities for two independent sources of thermal 
light with mean photon number eN. The expected visibilities for the thermal and 
Poissonian cases are compared in Fig. 12.51 These averages do not tend to one - as 
the visibility underestimates the degree of localisation when the prepared states are 
localised at two values of the relative phase, as discussed previously in Sec. 12.3.11 
However the general trend is clear. The cavity modes initially in thermal states tend, 
as in the Poissonian case, to a state which is perfectly correlated in relative phase 
while remaining unentangled. Although more photons must be detected to achieve 
the same degree of localisation when the initial states are thermal, the localisation 
proceeds very rapidly in both cases. 
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2.3.3 The consequences of non-ideal photodetection 

An important issue regarding the experimental feasibility of the interference proce- 
dure for preparing states with a well-defined relative phase is the problem of errors in 
the photodetection process. There are two types of errors — detector inefficiency (or 
photon loss) and dark counts (for a recent review of the current state of technology 
see |Migdall04| ) . In general, dark counts are more serious if they occur unpredictably 
in the experiment, and to incorporate them in calculations it is necessary to mix over 
the outcomes in the ideal case consistent with what is actually registered. However 
detector inefficiency presents the greater problem technologically (typically 20% at 
optical frequencies). An inefficiency detector can be modelled as in |Yuen83| . For 
the analysis of the interference procedure here where one of the detectors has effi- 
ciency T] < 1, the ideal detector assumed previously is replaced by a beam splitter 
with transmittivity and reflectivity a/1 — t], which couples the ancilla to an addi- 
tional vacuum mode, and with the primary output port now monitored by a perfect 
detector. The dynamical evolution during the procedure causes the additional mode 
to accumulate some population which is not measured, and the mode is traced over 
at the end. 

In fact photons lost from the apparatus due to detector inefficiency cannot affect 
the process of localisation directly. Instead, the beam splitter which combines the 
ancillae ensures that the origin of the lost photons from either cavity is uncertain, 
and the situation is equivalent to one where both cavities leak separately into exter- 
nal modes not involved in the measurement process. The localisation of the relative 
phase depends on those photons which are registered by the detectors in the same 
way as for the ideal case, while the photon loss reduces the total number of photons 
in the final state. For the case of large initial number states this situation has been 
studied numerically by M0lmer, as discussed in Sec. 12.1.41 The relevant simulation 
assumes loss rates out of the cavities equal to the leakage onto the beam splitter, 
a situation M0lmer identifies as equivalent to 50% efficiencies for both detectors. 
As another example, the calculations in Sec. 12.3.11 for the case of Poissonian ini- 
tial states can readily be extended to incorporate detector inefficiencies and, as is 
straightforward to see, the relative phase localisation conditioned on a particular 
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measurement outcome is as for the ideal case, although the probabilities for specific 
outcomes a fixed time after the start are altered. Finally it should be pointed out 
that postselection provides one potential aide when the photodetectors are subject 
to errors. For example, one can reject all runs of the experiment where the total 
number of detections is significantly different from what would be expected in the 
ideal case. 
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Advanced Topics on Localising Relative Op- 
tical Phase 



This chapter extends the programme of Chapter |21 on locahsing relative optical 
phase and the canonical interference procedure explained in Sec. 12.11 and presents 
original work on several topics. Sec. 13.11 formally proves the assertion that the 
action of photodetections during the interference procedure can be described by 
Kraus operators a ± b (where a and b are the annihiliation operators for the two 
modes). The full form of the measurement operators incorporates also the leakage 
from the cavity modes. The identification of the leakage parameter e with time is 
explained. 

Sec. 13.21 looks at the localisation of the relative phase for initial states with very 
different intensities at each mode, focusing on the example of initial Poissonian 
states. This example is treated along similar lines to the symmetric case discussed 
in Sec. 12.3.11 of Chapter El The evolution of the localising scalar function and the 
probability distribution for different measurement outcomes are found to be very 
different in the highly asymmetric case compared to the symmetric one. The issue of 
whether there are preferred values for the localised relative phase in the asymmetric 
case is discussed. Intuitively one might expect the localisation of the relative phase 
to be slower for increasingly asymmetric initial states, and this is verified here by 
looking at the evolution of the average visibility. In fact the definition of the visibility 
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given in Sec. 12.21 of Chapter |21 must be modified for this case. It is observed that 
the hmit of sharp locahsation becomes harder to attain for increasing asymmetric 
initial states. 

The transitive properties of the locahsation process are considered in Sec. 
Specifically the localisation process is assumed to act pairwise on three modes ini- 
tially in Poissonian states which can have different intensities. The conclusions are 
expected to hold generally. It is concluded that the localisation process acts inde- 
pendently of prior phase correlations with further systems, and that the localised 
quantum phases have the same transitive properties as classical phases. However 
the localisation process does affect phase correlations with external systems through 
(asymmetric) photon depletion. It is pointed out that for three modes with well 
localised relative phases loss of one of the systems does not disrupt the phase corre- 
lation between the remaining two. 

The possibility of using the canonical interference procedure for linear optical 
and classical feedforward based state engineering is suggested in Sec. 1114^ studying 
in detail the "addition" of photon number states. First the situation is analysed 
wherein two Fock states are combined at a beam splitter and one of the output 
ports is measured by a photodetector, to yield a larger Fock state at the free out- 
put port. An improved method of addition is then suggested, motivated by the 
fact that if the input states are phase locked, all the photons can be directed to 
one output port of the beam splitter using only an additional phase shifter. The 
probability for detecting the vacuum at the free out port, given input Fock states 
which are identical, can be doubled compared to the first method at the cost of one 
photon used to partially localise the relative phase, and almost tripled at the cost 
of two. As an aside it is pointed out that relational Schodinger cat states produced 
from large number states as in Sec. 12.1.31 of Chapter |2l can easily be converted to 
"NOON" states whenever the relative phases of the cat components are different 
by approximately vr. Hence the analysis in this section points towards a possible 
new route to generating NOON states based on procedures which establish relative 
phase correlations. 

Sec. 13.51 discusses the application of states prepared from initial Poissonian or 
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procedure 

thermal states as in Chapter |2l with well localised relative phase and a fixed total 
energy, for fundamental tests of superselection rules (as published in |Cable05j ). 
In algebraic quantum field theory the existence of absolute conservation laws, and 
associated superselection rules forbidding the creation of superpositions of states 
with different values of the conserved quantity, is taken to be true axiomatically. 
A less absolutist and more operational approach was initiated by Aharonov and 
Susskind jAharonov67j , who suggested that the forbidden superpositions can in fact 
be observed provided that the apparatus used by an observer are prepared in certain 
special states. The states suggested by Aharanov and Susskind are not particularly 
realistic. Here mixed states with well localised relative phase, which are much 
more experimentally feasible, are presented as alternatives which can reproduce the 
desired effects with no loss due to the lack of purity. To end the chapter. Sec. 13.61 
suggests possible future calculations following the programme set out in Chapters El 
and El 

3.1 Derivation of the measurement operators for 
the canonical interference procedure 

In the canonical interference procedure the effect of photodetections on the two 
cavity modes is given by Kraus operators Kl and proportional to a ± 6, where a 
and b are annihilation operators for the cavity fields, as assumed by fiat in Sec. 12.1.11 
of chapter O This result assumes ancillae which are initially in the vacuum state, 
and for simplicity the phase shifts in the apparatus are taken to be 0. In this 
section a mathematical derivation of these measurement operators is presented, and 
in addition the effect of leakage from the cavities is fully accounted for. The action 
of the interference procedure on an arbitrary tensor product of Glauber coherent 
states is evaluated, and corresponding operators for the full evolution of the system 
are deduced. In fact this is all that is required. The same measurement operators 
must also be valid for arbitrary initial states for the cavities, pure or mixed, since 
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coherent states form a complete basis and all the operators involved are linear.^ 

The canonical localisation process is described in detail in Sec. l2.1l of Chapter|21 
and is summarised in what follows. Each cavity mode couples to an ancilla mode, 
initially the vacuum, according to a linear mode coupling with small parameter e. 
The ancillae are combined at a 50 : 50 beam splitter, and the output channels are 
measured by photodetectors with r photons detected at the "right" detector and 1 
at the "left" detector. Mathematically the steps of this process are as follows.^ An 
initial product of coherent states |a)A® \P)b for the cavity modes and ancilla modes 
|0)a ® \0)b external to the cavities are first transformed as, 

|0),®|0),®|a)^®|/3)^ 
— ^ I v^a)^ ® I ® I Vl - ea)^ ® | Vl - ep)^ , 

due to the leaking of a small fraction of the cavity populations as determined by e. 
The 50 : 50 beam-splitter acts on the ancillae and they evolve to, 

K"^) ). H ^ (^) ). ^ ^ I ^"^^ 

Finally r photons are measured in the ancilla labeled a, and / photons are measured 
in the ancilla labelled b, corresponding to action of the operator \r){r\a ® \l){l\b, 
leading to the final state. 



Overall the measurement process involves decay of the cavity field amplitudes in 

^The general form of Kraus operators for a system which is coupled to an ancilla is derived in 
Section 8.2 of INielsenOO) . It is assumed there that the system and ancilla are initially uncorrelated 
with the ancilla in some state |eo) at the start. The system and ancilla then couple under some 
unitary process which denoted U. At the end the ancilla is projected onto one of the members of a 
complete orthonormal basis {|efe)} for the ancilla. The Kraus operators for this entire process are 
shown to be of the form (efc|J7|eo). The derivation in this section identifies the operators {ek\U\eo) 
relevant to the canonical interference procedure. Note that strictly speaking the discussion in 
[ NielsenOO] concerns a general system described by a finite dimensional Hilbert space, whereas in 
the current problem the state space of each mode is infinite dimensional. 

^It should be emphasized that this heuristic treatment of the dynamics is expected to be valid 
only when e, also the fraction of the cavity populations which has leaked into the ancillae, is small. 
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3. 1 Derivation of the measurement operators for the canonical interference 

procedure 

addition to the action of the operators a ± 6 on the cavity modes, and there are 
coefficients relating to normahsation. 

An obvious candidate for an operator causing the decay of the amphtude of a 
coherent state \a)A would take the form exp — (/t/2) iV^t , where Na is the number 
operator and k is a decay rate and t a time parameter. In a continuous measurement 
model of photon detection such an operator would generically correspond to a tem- 
poral evolution conditioned on no photon detection (for a review of quantum jump 
methods refer |Plenio98] and references therein). Expanding the coherent state \a)A 
in a basis of Fock states, it may be deduced that a suitable operator taking \a)A to 
\\/l — €a)A is of the form e^Ainv^T^_ (detail. 



For two modes with annihilation operators a and b the appropriate operator would be 
g(ata+btb) In v^T^^ Identifying a potential time parameter. In y/1 — e = — | (e + O (e^)), 
and hence e ~ nt, provided |e| ^ 1.'^ 

A complete measurement operator can now be proposed which has the correct 
action on |0)a ^ \0)b ^ \a)A ® \i3)b, 



1,(01 J0(o|.e.p {4 («W..) H. (1-4 ySp/^ . 

(3.1) 

Only the action of the procedure on the cavity fields is of interest so the initial 
component |?")(0|^ 1 0(0 1 6 Ki^iy be dropped. As for the additional coefficients: the 
^^/(l+Vj^/~(J\rT) is the root of the number of sequences of detections with / "left" 
counts and r "right" counts; the a/c'"^^/ {l+r)\ can be understood as the root of a 
factor relating to the probability of having a total of / + r detections. 

It is worth here emphasizing that this heuristic treatment of the dynamics leads 
to measurement operators with components, which can be identified as "jump oper- 
ators" and a "decay operator", which are also the principle elements of a quantum 
jump type analysis. However the latter analysis is based on many additional approx- 

^An identification of e and time t could be deduced just by considering dimensions. However 
the full expression i c>c — In y^l — e is not so obvious. 
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imations. Specifically the reservoir modes are assumed to be multimode (associated 
for example with a continuous spread of wavelengths and different directions of prop- 
agation). In addition they are assumed to be strictly Markovian — the reservoir 
correlation time is taken to be very much smaller than the typical relaxation time 
of the "principle system" . The interaction of the principle system with the reservoir 
is supposed "weak" over short times — so that all but the leading order interaction 
terms can be discarded (the Born approximation). One assumption which is shared 
with the heuristic analysis is that the reservoir and principle system are uncorrelated 
initially. Furthermore in both cases it is assumed that the governing Hamiltonians 
are consistent with the rotating wave approximation (for which components of the 
interaction Hamiltonian which do not explicitly conserve energy are discarded). In 
conclusion, the treatment of the dynamics in this section is seen to be reliable for 
short times despite its simplicity. 

3.2 Analysis of the canonical interference proce- 
dure for asymmetric initial states 

This section considers key features of the localisation of the relative optical phase 
between two cavity modes which undergo the canonical interference procedure when 
the initial intensity of one of the modes is very much greater than the other. Specif- 
ically the methods employed in Chapter |21 are extended for the case of initial asym- 
metric Poissonian states, an important case since — is as much as the state of 
laser light may be taken to be Poissonian — lasers with a wide range of intensi- 
ties is readily available in the laboratory. More generally this discussion is relevant 
to a wider set of issues — understanding the situation when a microscopic system 
is probed by a macroscopic apparatus; analysing reference systems composed of 
quantum resources; and explaining the emergence of classical behaviour in closed 
quantum systems. 

Suppose that two cavity modes undergo the canonical interference procedure 
described in Sec. 12.11 of Chapter |2l but now the cavities are populated initially 
by Poissonian states with different mean photon numbers N and M respectively. 
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Following the analysis in Sec. l2.H.n of ChapterEl the initial states of the two cavities 
may be expanded in terms of coherent states, 



pi = ^n„(iV)|n)(n|®^n„(M)|m)( 
1 



m\ 



n=0 m=Q 

ded<\)\ci)l^Q\ ® 1/5) 



47r2 



^0 



where a = Vn expiO and /? = Vm expicf), and n„ (/i) = /i^e^^'/n! is a Poissonian 
probability factor with mean /i. After / detections at the left detector and r at the 
right detector the state of the two cavities has evolved to, 

, ^ /■^Y'"^^^'/'^-erjv+M)e'+- fa+Py f a*+f3* y f -a+(3 ^' f -a*+(3* ^' 
Jo Jo 



27r27r l\r\ \ ^2 J \ ^2 J \ V2 J \ V2 



X \VT^a){Vl^a\0\Vl^p){VT^P\, (3.2) 

where e denotes the leakage parameter for both cavities and takes a value very much 
less than 1. p' here is subnormalised and the value of the trace is the probability 
for the specific measurement outcome at the two detectors. The scalar function 
\Ci^r (^)l^; where A = — 6*, is defined by. 



P = II \Ci,r {A)f \VT^ea){VT^ea\ ® |v^/3>(v^ 

Jq Jq 2tt 2tt 

and takes the form, 

\Ci,r{A)f oc (l + i^cosA)'■(l-i^cosA)^ (3.3) 



where R is defined as the ratio 2yNM / (iV + M) , which must take values between 
and 1. -R is 1 when N = M, 0.94 when N is twice M (or vica- versa), 0.57 
when the factor difference is 10, and 0.20 when the factor is 100. Noting that 
\Ci^r {A)f = \Cr,i (~^)|^; it IS seeu that the probability of I and r counts at the left 
and right detectors, and the degree of localisation of the relative phase parameter 
for the final state, are identical if the outcomes at the two detectors are swapped 
around. 
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Figure 3.1: The normalised scalar function |Ci_r.(A)p for a total of 15 photocounts and 
the number I of "left" detections going from to 7. Initial Poissonian states are assumed 
with the intensity in one mode 10 times the other (parameter R = 0.57). The dashed 
curves peaked at Aq = correspond to Z = 0, . . . , 3 and have progressively larger spreads. 
The solid curves are for I = 4, ... ,7 and are peaked progressively further from A = 0. The 
curves for the case of more left then right detections would be centered about A = vr and 
related symmetrically to the ones shown. 



For the highly asymmetric case when R takes smaller values, the \Ci^r (^)l^ are 
qualitatively different to the case of symmetric case R = 1 and are as follows. If 
the counts are predominately at the right detector, |Cz,r(A)p is peaked around the 
value Ao = for the relative phase. Given a fixed total number of detections 
/ + r, the localising function is progressively less well localised as the number of left 
counts increases, leading to a qualitative change when / and r are close in value. 
Specifically, for the cases when r—R{l+r) < I < r, |C/_r.(A)p is peaked at two values, 
Aq = ± arccos ^^^^ ■ Differently here from the cases of initial symmetric Fock or 
Poissonian states (discussed in Sec. 12. 1.31 and Sec. l2.3.I] of Chapter 12)), |Q^r(A)p is 
non-zero at A = and a partition of the final state into two components localised 
at +Ao and — Aq is harder to justify. When the I > r the pattern of the |C/,r(A)p is 
reversed, and when the detections are predominately at the left detector |C/^,.(A)p is 
peaked at Aq = vr. The situation is illustrated in Fig. 13. 11 for the case that r > I, and 
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3.2 Analysis of the canonical interference procedure for asymmetric initial states 

a total of 15 detections at both detectors. The initial state is chosen with R = 0.57 
corresponding to the case that one cavity begins with 10 times the intensity of the 
other. Dashed curves are for I = 0, ... ,3 and solid ones are for I = 4, ... ,7. 

For a fraction of approximately (1-R) of the measurement outcomes, the relative 
phase is localised at or vr, and the question arises as to whether these are "preferred 
values" when the initial states are highly asymmetric. To answer such questions it 
is necessary to investigate the behaviour of Pi^r, the probability distribution for 
different detection outcomes, for small values of the parameter R. A full expression 
for the probability Pi ^ of I "left" and r "right" detections is given by. 



tr ip') 



n,.,(eiV + eM)| (^) (-1)' 




(3,4) 



where (*) abbreviates the combinatorial factor (s + t)! / (s!t!), and p^egendre (jg^otes 
the Legendre polynomial of degree i given by. 

The Poissonian factor Hi^r {^N + eM ) depends on the leakage parameter and the 
sum of the initial cavity intensities, and controls the probability for the total number 
of detections l+r. Evaluating P/^^ numerically reveals the following. The distribution 
changes qualitatively as R is reduced in value. For small values of R, The situation 
is essentially reversed from the symmetric case R = 1 (refer Fig. 12.31 of Chapter 12)). 
Given a fixed total number of detections, outcomes with detections at one detector 
are found to be least likely. Instead, the probabilities increase to a maximum at 
I = r, (or this peak is split and the most likely events are at two points for which 
r ^ /). 

Turning attention now to the total probability of localisation of the relative 
phase at Aq = 0, vr, let Po,7r(-R, -D) denote the total probability of such events for 
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a fixed total number of detections D = I + r, and wliere tlie tfie Poissonian factor 
Ui+ri^N + eM) is divided out. Events witli localisation at 0,7r are determined by 
the inequality, 

\l-r\> R{l + r). (3.5) 

For the case of symmetric initial states, -Po,7r(l, D) = 2r(D + 0.5)T{0.5)/'kT{D + 1), 
which is well approximated by 2/ \fnD. Evaluating Pq j^(R, D) numerically for other 
values of R reveals the following. For cases when the initial state is asymmetric but 
R is not too small, Po^n{R,D) follows the same general trend as in the symmetric 
case although it is weakly oscillatory and can be bigger or smaller than Po,7r(l,-D). 
In fact, even for R as small 0.57, when the intensities of the cavity modes are 
initially different by a factor of 10, Po,7r(0.57, D) and Po,7r(l;-D) take similar values 
for D up to 100. However, for more extreme values of R, Po^n{R, D) is found to be 
consistently greater than Po,7r(l, -D). For example, for R = 0.2, Pq^t,{0.2, D) is very 
roughly double the probability for symmetric case for D up to 50. Going to still 
smaller values of R, Pq .„{R, D)I Pq .„{\, D) is found to increase both with increasing 
detections D and R tending to 0. In conclusion, in the limit of highly asymmetric 
initial states, the increasing probabilities for events with / ^ r is less significant 
compared to the increase in the proportion of measurement outcomes satisfying 
relation ()3.5|) . 

Thinking now about applications to understanding the transition between quan- 
tum and classical behaviour in the natural world (or even decoherence in controlled 
laboratory experiments), the canonical interference procedure could be viewed as 
a simple prototype for a natural process involving essential elements of indistin- 
guishability (provided by coherent combination under a linear mode coupling) and 
measurement. It is of interest then to consider the localisation of the relative phase 
after mixing over all possible measurement outcomes. For the case of highly asym- 
metric Poissonian initial states it might be expected on the basis of the discussion 
above that the averaged mixed state describes localisation of the relative phase at 
and TT. In fact this is not so, which can easily be demonstrated for the general 
case of initial states of the form Pa = j ^Pa{oi) and Pb = J ^Pb{I3) \l3){l3\. 

Here the modes are labelled A and P, |a) and \[3) are coherent states and PA(a) 
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and Pb{i3) are the Glauber-Sudarshan P representation |Glauber63| ISudarshan68j ) 
of the states. Under the action of leakage with parameter e into the external modes 
labelled Ext A and ExtB, 



Pi = PA®PB 



(g) I v^I^a )( v^T^a I ^ ® I )( 1 



B ■ 

Letting Ucomb denote the combined unitary action of fixed phase shifts and a linear 
mode coupling of arbitrary strength and duration on the external modes, and sup- 
posing a subsequent measurement of / and r photons in modes Ext A and ExtB, 
the final state of modes A and B is, 



P^^r = h^^PA{a)PBi/3) {l\{r\U,omi>W'a)V~^/3) 



Ext 



X \VT^ea){VT^ea\^0\VT^ep){VT^eP\^, 

where pi^r is subnormalised and tr (p;,r) is the associated probability. If the record 
of the photodetections is unavailable, the final state of modes A and B is given by 
mixing over all measurement outcomes and. 



oo 



Pf = H'-M^^ (3.6) 



Z,r=0 



Z,r=0 

However summing over the localising factors, 

oo 

^ |(/|(r| f/comb |Vea)|ye/3)E^J = 1- 



;,r=o 



Hence this complete mixing washes out the measurement-induced localisation. This 
result is in contrast to the situation when a state such as a Poissonian state is beat 
against the vacuum at a beam splitter. In this case a relative phase is acquired 
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depending on fixed physical properties of the beam sphtter, and something similar 
would happen for highly asymmetric states for which both modes are populated. A 
different point of view on the mixing process Eq. ()3.6|) is provided by the theory of 
decoherence as proposed by Zurek et al |Zurek03j . and is discussed in the Outlook, 
Chapter IHl 

Returning to the example of initial Poissonian states, another question concerns 
the speed of localisation for the symmetric versus asymmetric case. This issue is 
complicated by the fact that even a perfectly phase locked state. 



Plocked = j^^^ |^e^^)(v^e^^| ® I v^e^(^+^°))(v^e^(^+^«)| , (3.7) 



(where the basis states are coherent states) has a visibility less than 1 whenever 
the amplitudes \fN and \fM are different from each other (for the definition of 
the visibility refer Sec. 12.21 of Chapter El), and will lead to an imperfect pattern 
of interference. In fact, the visibility for piocked can easily be shown to be 
equal to the parameter R assigned above to products of Poissonian states with 
corresponding amplitudes. However, a combination of a single mode phase shifter 
and a beam splitter with reflection and transmission coefficients determined by the 
ratio ^JN/M is sufficient to arrange deterministically for one or other of the modes 
of Plocked to be in the vacuum state, provided the phase parameter Aq is known 
accurately. Hence, it seems reasonable to attribute perfect phase localisation to 
Plocked and to define a modified visibility taking values between and 1. In what 
follows, visibilities are computed for the partially localised states Eq. ()3.2p generated 
from initial products of Poissonian states by the canonical interference procedure, 
and are adjusted by a factor 1/R. It should be noted that in this example the 
leakage and measurement processes do not alter the relative intensities of the two 
modes, and that in the limit of sharp localisation, the final state will be of the form 
of 

plocked- 

The symmetry of the localising scalar functions |Ci,r(6', 0)P, as illustrated in 
Fig. 13.11 allows a simple formal expression for the visibilities to be written down 
in terms of the probability distribution P/ for corresponding detection outcomes. 
When / = r the visibility can be shown to be 0. In the case of more right than left 
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(a) Evolution of the expected visilibility for 
increasingly asymmetric initial states 
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(b) Evolution of the expected visibility 
for events localising at orTi 
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Figure 3.2: Evolution of tlie expected visibility ^ Pi,r^i,r with increasing values of 
e{N + M) for initial asymmetric Poissonian states with intensities N and M in each mode. 
Curves are for R = 0.94, 0.57 and 0.2 and increase more slowly for smaller values of R. 
In (a) all possible detections outcomes are included in the averaged visibility. In (b) 
only outcomes with localisation at a single value of the relative phase are included in the 
average. 



detections, r > I, the intensity /(r) in the definition of the visibility (refer Sec. 12.21 
of Chapter 12)) is maximised and minimised for phase shifts of r = 0, vr for which it is 
proportional to (r + l)Pi^r+i and (/ + l)Pi+i^r respectively. The situation when / > r 
is analogous. The following expression for the visibility, including the additional 
rescaling factor l/R, 



(r + l)P;,,,+i -(/ + !) P/+1,. 



l+l,r 



is valid for any outcome at the detectors. 

The evolution of the expected visibility J2i r Pi,r^i,r as a function of e(iV + M) 
is plotted in Fig. 13.21 (a) for the cases R = 0.94, 0.57 and 0.2. The average rate 
of localisation of the relative parameter is seen to be progressively slower for more 
asymmetric initial states. For the slowest case R = 0.2 the expected visibility attains 
85% the value for the fastest case R = 0.94 by e{N + M) = 100, while the cases 
of i? = 0.57 and R = 0.94 converge over the range shown. More analysis would be 
required to clarify whether the localisation can be as sharp (as in the symmetric 
case) for small values of R given enough detections. However, the slower rate of 
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localisation for highly asymmetric initial states would in any case necessitate a large 
total intensity for the initial states if the limit is to be attained. Another feature of 
Fig. 13.21 (a) is that the expected visibilities increase to a value very much less than 
one — in fact the maximum is 0.63 on the range shown. This reflects localisation 
at two values of the relative phase. In Fig. 13.21 (b) the averaging for the visibility is 
restricted to outcomes for which there is localisation at single values, which can be 
or vr (in other words outcomes satisfying Eq. ()3.5|) ). and account is taken of the 
total probability of such events. The expected visibilities for these events is 0.98, 
0.94 and 0.81 by e{N + M) = 100 for R = 0.94, 0.57 and 0.2 respectively. 



3.3 Transitivity of the canonical interference pro- 
cedure 

This section considers the transitive properties of the process of localisation under- 
lying the canonical interference procedure, introduced in Sec. l2.1l of ChapterEl The 
following examples assume that the states of the optical modes prior to the localisa- 
tion processes are Poissonian, which can readily be prepared in the laboratory, but 
the conclusions are expected to apply generally. 

Suppose that two systems, modes 1 and 2, have been prepared in a phase locked 
state by an interference experiment in which several photons were detected, and are 
brought together with a third system, mode 3, in a Poissonian state with no phase 
correlation with either mode 1 or mode 2. The state of the first two modes is, 



(1.2) 



2^ d9 
2^ 



and of the third mode is. 



ie 



(3.9) 



(3) 

Pi 



27r 



2n 



me'- 



i<j> 



(3.10) 



where the basis states are coherent states with amplitudes ^/iVj and j = 1,2,3 
is the mode label. Suppose now that the interference procedure is performed on 
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modes 2 and 3, with / and r photons detected at the "left" and "right" photocoun- 
ters, as described in Sec. 12.11 The state of the three systems is transformed as 
pf''^^ ® pf^ —>■ p' as follows, 

p' oc (a2 + asY (a2 - 03)' p^^'^^ (g) pf^ (^al - al^ (^al, + 4 j 

r2iT r2iT 



Jo 



dOdcj) 
27r27r 



Ci,r (^,0)r 



j(9+Ao) 



1 



i<j) 



, (3.11) 



where p' here is assumed to be subnormalised so that tr(p') gives the probability of 
the specific outcome at the detectors. The scalar function jCj^rP encoding informa- 
tion about the localisation of the relative phases takes the form, 



Lr I 



OC 



2r 



21 



(3.12) 



where A = (p — 6 — Aq. For example if |Ci^,.p is sharply peaked around the value Ai 
for A, the final state of the three systems is. 



p oc 



2lT 



de 
2^ 



]V2e'(e+Ao)^^y]^gi(e+Ao) 



/]^gi(e+Ao+Ai) y ^y]^gi(e+Ao+Ai) 



The relative phase of mode 3 relative to mode 1 has evolved to Aq + Ai. If the 
second system is lost — i.e. mode 2 is traced out — the final state of modes 1 and 
3 is, 



2n 



de 
2^ 



i(6l+Ao+Ai) 



and the phase correlation remains unaffected by the loss. 

The previous example demonstrates that the localised relative quantum phases 
have the same additive properties as classical phases between different systems. 
Furthermore, the interference procedure on the second and third systems is inde- 
pendent of the prior phase correlation with the first system. This is clear from the 
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dependence of jC^^rP on A = — — Aq and the fact that, 

''^P^ = L Jo ^^^^''^"^-^-^'^^ =1 Jo • 

(3.13) 

tr(p'), determining the probabihty for different measurement outcomes of the in- 
terference procedure, is the same here as for an initial product of Poissonian states 
with corresponding mean photon numbers A^2 and N3 uncorrelated with a further 
system. 

In practice however estabhshing a relative phase between two systems will have 
some effect on phase correlations with external systems due to the effects of photon 
depletion. With reference to the previous example, if the amplitudes of modes 2 
and 3 are reduced by a factor Vl — e during the interference procedure, the relative 
amplitude of modes 1 and 2 changes as ^/Ni/N2 ^ (l/v^T^) This 
does not affect the value Aq of the relative phase between modes 1 and 2, but does 
change, for example, the behaviour of the two modes if combined on a beam splitter. 

The independence of the canonical interference procedure of phase correlations 
with external systems makes it interesting to consider what happens in a three mode 
system where, as previously, mode 3 is initially Poissonian, and an interference 
procedure is performed on modes 2 and 3, but where differently, modes 1 and 2 
begin in a state where the relative phase is localised at two values iAg. Such states 
are produced in an ideal interference experiment (refer Sec. 12.1.31 of Chapter |2]). It 
is easy to see that after the second interference procedure on systems 2 and 3, which 
for simplicity is assumed to determine the value for the relative phase, the final 
state of the three mode system is. 




Interestingly if the second system is lost (i.e. if mode 2 is traced over), modes 1 and 
3 are left in a state localised at two values, but if the first system is lost (i.e. if mode 
1 is traced over), modes 2 and 3 are left localised at one value of the relative phase 
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0. It should also be pointed out that for a well locahsed state / |^|a;)(a;| \P){P\, 
where a = a/(1 — e)A''ie*^ and /5 = a/(1 — e)N2e^^^^^°\ the expected numbers of 
photons in each mode — (1 — e)A^i and {1 — e)N2 — are independent of Aq, and 
hence the effects of photon depletion in an interference experiment involving one 
mode of pair of systems in a state having two components with different localised 
relative phases will in general be the same on both components. 

3.4 Application to engineering large optical Fock 
states 

This section presents an example of how the canonical interference procedure in- 
troduced in Sec. 12.11 of Chapter El can be utilised in schemes for quantum state 
engineering based on linear optics and classical feedforward of measurement re- 
sults. The example considered is the generation of large Fock states from a source 
of single photons by a probabilistic process of "addition". Large photon number 
states are of interest for several reasons. They have applications to interferometry 
with phase uncertainty reduced to the Heisenberg limit, considered theoretically in 
|Holland93l IKim98j . Second, the coherence properties of Fock states are of funda- 
mental interest, particularly in relation to those of Glauber coherent states. Finally 
relational Schrodinger cat states, generated by applying the canonical interference 
procedure to pairs of Fock states, can easily be converted to "High NOON states" 
— states of the form |0)2 + 10)^^ |A^)2) in a basis of Fock states where is 

assumed to be large — whenever the two components of the cat have relative phases 
approximately vr apart. This conversion is explained at the end of the section. For a 
recent discussion of some important applications of NOON states, in particular for 
interferometry beyond the shot-noise limit, see |Leen5j . 

The problem considered here is how to create large Fock states given (ide- 
alised) linear optical components — beam splitters, photon detectors, phase shifters 
and a highly efficient source of single photons (for recent progress on sources see 
|Grangier04| and on detectors see |Rosenberg05| ). One obvious and simple but in- 
efficient method is as follows, suggested by the well-known Hong, Ou and Mandel 
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dip experiment |Hong87| , whereby two uncorrelated and identical photons, simul- 
taneously incident on the input ports of a 50 : 50 beam splitter, must both be 
registered at the same output port. Starting with a product of one photon states 
coherent combination of the optical modes at a 50 : 50 beam splitter,^ yields 
a superposition of states with two photons in a single mode, 

y|{-|2)j0)2 + |0)j2)2}. 

Assuming the second output port is monitored by a (non-discriminating) photode- 
tector, a two photon Fock state is created at the first output port with probability 
1/2. This process could be iterated to yield successively larger Fock states. For 
example if two two-photon states are combined on the beam splitter, 

|2)i 12)2 10)2 - 2 |2), |2), + v^|0), 14)2 }, 

in a basis of Fock states, and if the second mode is monitored by a photon-number 
discriminating detector, the output state in the first mode is |4) with probability 
3/8, and |2) with probability 1/4. For the case of an input state with the same 
number of photons in each mode it is seen that in half, or a little under half, of 
such operations the output is worse than the inputs. 

Looking in more detail it is helpful to determine a general expression for the 
action of a 50 : 50 beam-splitter on the state |A^)i|A^)2 with photons in each 
mode, which can be done as follows. Adopting a representation in terms of coherent 
states, 

where IItv (A^) denotes a Poissonian factor. Under the action of the beam splitter 

■^Fixing the phase convention for calculations in this section, the beam splitter transformation is 
taken to be exp (f — f 0^6), a and b denoting the annihilation operators for the modes labelled 
1 and 2. 
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this is transformed to 
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Expanding the coherent states in the Fock basis and making the replacement 

/oXf = /o7o"4¥)4^) [hfi-] (where ^ and ^ are the vari- 
ables of integration), vahd because of the 27r-periodicity in 6 and 0, 
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where A = (p — 6. The integral over the average of the phase variables enforces 
conservation of photon number A + B = 2N. As for the second integral the following 
result is required, 
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for an integral over the first quadrant, where r(-) denotes the usual gamma function. 
The integral over the phase difference variable A/2 over all four quadrants is seen to 
be four times this when A and B are both even, and otherwise since positive and 
negative valued contributions cancel. Overall {ip-^) is seen to be a sum over contribu- 
tions with even photon numbers for both modes. The final expression simplifies to 
the following (using the standard identity T (M + 0.5) = {2My.y/n/ {2'^^'^M\) vahd 
for M > 0), 
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(3.14) 



(in agreement with |Campos89| ) . When a photon number discriminating detector 
monitors the second mode the most probable outcomes are measurements of or 
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2N photons, and the least hkely measurement is N/2 photons (or (A^ ± 1) /2 pho- 
tons if odd). If a non-discriminating detector monitors the second mode only 
a measurement of photons would be useful, which happens with probability Pq 
where, 



Po = \{2N\,{0\,\^pf)\ 

{2N)\ 
~ 22F/Y12 ' 



(3.15) 



which well approximated by the expression 1 / V nN ( obtained using Stirling's for- 
mula for large factorial numbers). 

It is interesting to compare this result with the action of a lossless beam splitter 
on a product of coherent states |a)i|/?)2 or a phase locked state 
/o" ^e-'^^'^ 1/5)2, where a = v^e^^ and (3 = v^e^(^+^°) — products of co- 
herent states evolve deterministically to products of coherent states. A combination 
of a single mode phase shifter and a beam splitter with reflection and transmis- 
sion coefficients determined by the ratio \/Ni/N2 of the amplitudes of the coherent 
states suffices to get mode 2 into the vacuum state. This motivates an approach to 
"adding" Fock states based on the interference procedure discussed in Sec. 12.11 of 
Chapter 121 which localises the relative phase degree of freedom. Since the localisation 
in such a procedure is rapid, it would be expected that a substantial improvement 
in the success probability of the previous method is possible by sacrificing a small 
number of photons, and proportionally less for larger input states. In fact when the 
input states are one-photon Fock states there can be no benefit — a single photon 
Fock state could be obtained with certainty but is the same as the inputs. 

In detail a revised method of "adding" a product of Fock states |A^)i|A^)2 is as 
follows. Refer back to Sec. 12.11 of Chapter |21 for an explanation of the canonical 
interference procedure; the notation and assumptions here are the same. It is as- 
sumed now that: the canonical interference procedure is applied first and stopped 
after a fixed total number of detections W have been recorded; that the appara- 
tus is ideal i.e. phase stable; and that phase shifts throughout the apparatus are 
fully characterised so that full information is available about the value Aq of the 
(partially) localised relative phase. If the process is performed slowly then only 
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non-disciminating photocounters would be needed. If > 2 it is assumed that the 
process is interrupted half way through and that a phase shifter in the second arm 
(refer Fig. l2.1l of Sec. 12. 1|) is adjusted by 7r/2. This enforces locahsation at one value 
of the relative phase. Without this step a relational Schrodinger cat state would 
be obtained with high probability and it would not be possible to arrange for such 
large variation in the intensity at one of the output ports of a beam splitter applied 
subsequently as can be achieved otherwise. The state of the optical modes in a basis 
of Glauber coherent states is then, 
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where A = (f) — 6, N' = {1 — e) N (where e is a leakage parameter) and the scalar 
function C'°(^) is peaked at A = Aq and is normalised such that ('?/''°^|'?/''°'^) = 1. 

The next steps in the procedure are a phase shift of the second mode by — Aq + vr 
and coherent combination of the two modes at a 50 : 50 beam splitter, so that the 
final state is, 
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If a detector monitoring the second output port registers photons, 
onto the Fock state \2N — W)-^ as follows. 



is projected 
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where II2N-W [2 (1 — e) A^] is a Poissonian factor with mean 2 (1 — e) A^, and the 
variable of integration A/2 runs between and 27c. If a photon number discrimi- 
nating detector is available for monitoring the second output port then the output 
at the first output port need not be discarded when the detector fails to measure a 
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vacuum state; in general the yield of Fock states with more photons than the input 
states is expected to be considerably improved compared to the earlier method. 

Letting denote the normalised state of the two optical modes after the 

localisation procedure has been performed and after a phase shift of — Aq+tt has been 
applied to the second mode, \ipP^°'^) = exp [i (— Aq + tt) \tp'-°'^) where b denotes 
the annihilation operator on mode 2, and letting U (50 : 50) denote the beam splitter 
transformation, the amplitude for detecting the vacuum state is given in a basis of 
orthogonal Fock states by, 

(2Ar-iy|^(0|2f/(50:50)|^/'^"'°'^) 

(2iV-iy)/2 I oo f2N-w\^^'^ 1 
= (i) E g (_i)<--»--) (,|, (2iV -W-s\A 1^-) . 

(3.17) 

Consider the case when W = 1 and the interference procedure involves just one 
detection. If the photon is detected at mode 1 then clearly no additional phase shift 
is required. After normalisation, 

^^proc^ ^ \N-l) \N)-\N)\N~1) ^ 

Pq, denoting the success probability |(2A^ — l|i(0|2f/(50 : 50) evaluates to 

2 (2A^)!/2^^A^!^, exactly double the value for the first procedure (refer Eq. ()3.15|l ). 
and consistent with the Hong, Ou and Mandel phenomenon when = 1. If the 
localising detection occurs mode 1, an additional phase adjustment of vr is required, 
and Pq is the same. 

Consider now the case when W = 2, and the detections are both at mode 
2. The protocol demands a fixed phase shift of 7r/2 between two detections so 
that 1?/''°'^) oc (a + 6) (a + ib) |A^)-^ |A^)2- Other outcomes at the detectors work out 
analogously. After normalising, 

\r'^) = \f^\N-2\ |Ar)2 + il±a^ lAT-l), \N-l),+^.^^\N), \N-2),, 
1^ / V 4Ar-2 ' ^JAN-2 ^ ^ V 4A^-2 ' ' 
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peaked at a relative phase of Aq = —tt/A. Hence a phase shift on mode 2 of 57r/4 is 
required to maximise the success probabihty Pq, after which the state of the optical 
modes is, 
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Then the success probability Pq is given by, 
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A simpler expression valid asymptotically is, 
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having a fractional error of less than 5% for > 15, but it is invalid for small 
values of A^. Taking the full expression for Pq, the ratio of the success probability 
here with that of the method without the interference procedure Eq. ()H.15j) is 2.6 for 
the uninteresting case N = 2 (when the input and output states are the same size), 
rising to 2.9 for larger values for A^. These preliminary results for W = 1,2 suggest 
that sacrificing even more photons to achieve sharper localisation of the relative 
phase may be helpful for optimising the addition process for larger Fock states. 

One potential application of engineering large photon number states is generating 
relational Schrodinger cat states, using again the canonical interference process set 
out in Sec. 12. II of Chapter |21 and from these NOON states. This could be done as 
follows, continuing with the same notation and assumptions. Consider the case after 
several detections (a total number D) for an initial product of Fock states |A^)i|A^)2 
when the counts at the "left" and "right" detectors are equal. The final state of 
the two optical modes is a relational Schrodinger cat state with components having 
sharply defined relative phases of ±vr/2, and is of the form. 
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where I7) denotes a coherent state, 7 = I7I e*^ and I7I = y^N — (D/2). Apphcation 
of a phase shift of 7r/2 to the second mode and combination of the two modes at a 
50 : 50 beam sphtter yields a NOON state. Under the phase shift the relational cat 
state is transformed as, 

and the beam splitter transforms the state to, 

d^,-A.?i2eH^m |o)^ 127), + ^e-l-I^P^-('^/^)l I27), |0), 
^ |0), \2N-D), + exp (i|7P7r) \2N-D), |0), ^ ^^^^^^ 

where the last expression is a superposition of states with 2N — D photons in a 
single mode as required. In practice it would be necessary to limit the number of 
detections in the interference procedure as much as possible, leading to components 
with visibility significantly less than 1, and use the states generated when the counts 
at the left and right detectors are similar but not identical (when the difference of 
the localised relative phases is different from tt) so as to increase the yield. The 
resulting states would be expected to have a NOON state as a large amphtude 
component. Future work could clarify the usefulness of such approximate NOON 
states and optimise the suggested procedure. 

3.5 Relative optical phases and tests of superse- 
lection rules 

A conservation law makes operational sense (or nonsense) only when related to the 
procedures whereby the conserved physical quantities are measured. In particular, 
the frame of reference against which the measurements are made plays a crucial 
role. Certain frames of reference (e.g. position, atom number) are more in accord 
with everyday experience than others (e.g. "charge phase", "isospin phase"). This 
has perhaps more to do with the ground state of the universe (the electromagnetic 
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vacuum in particular) which acts as a readily accessible reference frame, than with 
any fundamental physical restrictions. 

A belief in absolute conservation laws leads to a belief in absolute superselec- 
tion rules. For example, in algebraic quantum field theory the existence of absolute 
conservation laws and associated superselection rules, forbidding the creation of su- 
perpositions of states with different values of the conserved quantities, is taken to 
be true axiomatically |Wick52t Haag96| (see also the discussion in section 2.7 of 



|Weinberg95| ). To illustrate how the more relational approach works consider quan- 
tum optics under the "rotating wave approximation" (refer for example |Gerry04] ) , 
taken to be equivalent to a strict superselection rule for energy under which the 
energy is additively conserved, and in the absence of any absolute phase refer- 
ence. Under these assumptions superpositions of states of different photon num- 
ber (such as Glauber coherent states of light) are forbidden, as are superpositions 
of non-degenerate atomic states {\g), |e)). Interaction Hamiltonians are explicitly 
excitation-conserving so, for example, the familiar Jaynes-Cummings Hamiltonian 
H = \e){g\ ® 0- + \g){e\ is allowed (where a and are annihilation and creation 
operators for an optical mode). For a discussion of Jaynes-Cummings dynamics see 
|Gerry04| . 

If asked how to operationally create and verify the existence of a superposition 
of atomic states of the form \g) + |e), a simple response would be to send the atom, 
initially in the ground state \g), through a cavity containing light in a large amplitude 
coherent state^ for an appropriate length of time, where the interaction Hamiltonian 
is H above. Measurement of the atom after exiting the cavity yields it in the ground 
state half the time. However this could be due to the atom being in a mixed state. 
To verify that a coherent superposition has been obtained the atom could be sent 
through a second cavity, also in a optical coherent state with a large amplitude, 
after which it would be found that the atom is always in the excited state |e) for 
example. This demonstrates that the atom was indeed in the coherent superposition 
\g) + \e) between the cavities. Note that this notion of a superposition assumes that 
the phase reference provided by the state of the first cavity is maintained. 
^Sufficiently large to drive Rabi oscillations between \g) and |e). 
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What Aharonov and Susskind (AS) note in |Aharonov67j is that the two cavities 
need not, in fact, be prepared in coherent states. A state of the form Eq. ()2.10|] 
discussed in Sec. 12.11 of Chapter |21 is operationally just as good as initial coherent 
states for the purposes of demonstrating a coherent superposition of atomic states 
using the procedure described above. ^ As noted previously, this state has fixed total 
energy, and thus there would be no violation of the conservation law globally. An 
objection to this test having physical relevance can be made along the lines that 
states of the form Eq. ()2.1()|1 are not easy to come by in nature, and in fact are 
highly entangled. However the results of Sec. 12.3.11 and Sec. 12.3.21 of Chapter |21 
demonstrate that a mixed state of the form Eq. ()2.14j) would do just as well for the 
operational demonstration of coherent superposition envisaged by AS. Such mixed 
states are much more easily preparable, and would seem to conform more closely 
with the type of reference frame states that observers typically prepare. 

There are some interesting ways this topic could be extended in future work. 
One assumption made by AS in their calculations is that the cavity modes are 
prepared with large intensities. It would be of interest to check how effective a 
phase locked state with only a small total number of photons would be. Also it would 
be useful to check what happens when the interaction Hamiltonian is taken to be 
H = {a + a''){\g){e\ + \e){g\), which is derived from the electric-dipole approximation 
without also making the rotating wave approximation. Perhaps the easiest way to 
create the phase locked state Eq. ()2.10|) is by splitting the light from a single laser 
at a beam splitter, rather than by phase locking two independent light sources. 
However the process suggested by AS to demonstrate coherent superposition could 
also be used to investigate features of the (partially) localised relative optical phase 
induced in states prepared using the canonical interference procedure of Sec. 12.11 of 
Chapter El 

®AS actually considered the creation of a superposition of a proton (equivalent to \g)) and a 
neutron (equivalent to |e)) using coherent states — or otherwise — of negatively charged mesons 
(equivalent to photons), under a Hamiltonian identical to the Jaynes-Cummings one. 
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3.6 Possible extensions 

It is clear that the programme of investigation in Chapters El and El can be extended 
to address further questions about locahsing relative optical phases. For example, 
with reference to the canonical interference procedure, further calculations could 
be done for a setup for which the ancilla modes are combined at an unbalanced 
beam splitter, for which the reflection and transmission coefficients have different 
magnitudes. The indistinguishability that the beam splitter provides is central to 
the process of localisation, and it might be expected here that the localisation is 
slower, or only partial even after very many detections. It would be interesting to 
clarify how the process localisation differs for the two cases of asymmetric initial 
states (with different intensities in each mode) together with a 50 : 50 beam splitter, 
and symmetric initial states together with an unbalanced beam splitter. Another 
question asks what happens when there is uncertainty in the path lengths of the 
apparatus during the interference procedure. In this case the Kraus operators cor- 
responding to photon detection take the modified form d ± [exp (ix)] b (where a and 
b are the annihilation operators for the two modes, and x is a position operator for 
the path length), and the process of localisation acts on the sum of the path length 
and the relative optical phase. 
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4 



Interfering Independently Prepared Bose - 
Einstein Condensates and Localisation of the 
Relative Atomic Phase 



This chapter looks at the interference of two Bose-Einstein condensates (EEC's) 
which have been prepared separately and therefore have no prior phase correlation, 
a process in which localisation of the relative atomic phase plays an important 
role. The focus is on spatial interference experiments in which the condensates 
are released from their traps and are imaged while falling, as they expand and 
overlap. High contrast interference patterns have been observed experimentally 
|Andrews97[ IHadzibabicn4] . Javanainen and Yoo were the first to provide a detailed 
numerical demonstration that such interference can be predicted without assuming 
that each condensate has an a priori phase |Javanainen96bj . Several studies have 
been done by other authors, and |Cirac96l INar aschewskiQ 6 1 IYoo97t IHorak99| should 
be highlighted, though little was done analytically (an exception being |('astin97] ). 

This chapter presents and extends work published in |Cable05| . investigating 
the key features of the localisation of the relative atomic phase as a pattern of 
spatial interference is built up by a series of atomic detections. The analysis has 
close connections with the study of localising relative optical phase in Chapter |21 
Most surprisingly it is shown that the localisation of the relative atomic phase takes 
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place on the same rapid time scale as in the optical example, and much faster 
than the emergence of a well-defined pattern of spatial interference in the numerical 
simulations of Javanainen and Yoo and others. 

The chapter proceeds as follows. The conceptual difficulties associated with 
assigning a BEC a coherent state macroscopic wavefunction with a definite a priori 
phase are reviewed in Sec. 14.11 The spatial interference of independently prepared 
EEC's is analysed in depth in Sec. 14.21 A simple measurement model is presented in 
Sec. 14.2. 11 based on the approach adopted in |.Tavanainen96bj . The visibility defined 
in Sec. 12.21 of Chapter |2l which can be used to quantity the degree of localisation 
of the relative phase, is found to have a simple interpretation here in terms of the 
probability density for detection of an atom at a particular position. The case 
of initial Poissonian states is treated analytically in Sec. I4.2.2[ and the analytical 
predictions are compared with simulated results in Sec. 14.2.31 Sec. 14.2.41 discusses 
the issue of what information is available in principle about the localisation of the 
relative atomic phase from a given pattern of spatial interference. Finally in Sec. 14.31 
a summary is given of a recent experiment and several different theoretical proposals 
concerning methods which probe the relative atomic phase by optical means, thereby 
performing a non-destructive measurement. 

4.1 The standard story and conservation of atom 
number 

A common and useful description of EEC's makes use of a coherent state macro- 
scopic wavefunction for the condensates. Such a description is generally justified by 
invoking standard stories about symmetry breaking (see for example Chapter 2 of 
|Pitaevskiin8j and |,Tavanainen86j ) . However there are several reasons to be suspi- 
cious of this common treatment. The first is that it requires a description of the 
state of a EEC in terms of a coherent superposition of states with different atom 
numbers. If atom number is to be conserved then such a description is tricky to 
justify. It is often argued that the EEC is surrounded by a thermal cloud with which 
it is exchanging atoms, and thus the atom number is undetermined. However such 
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a process leads only to a mixed state for the BEC and not a pure coherent state. 
Secondly, the symmetry breaking is generally invoked by the addition of auxiliary 
fields with no clear physical relevance (see for example jNar aschewskiQ 6 j for a de- 
scription and criticism of an approach based on Bogoliubov auxiliary fields). An 
assumption of an a priori symmetry breaking also implies the philosophically pre- 
carious position of writing down quantum mechanical states containing in principle 
unknowable parameters - the absolute atomic phases^. Finally the most striking 
demonstrations of coherence for EEC's come from interference experiments such 
as reported in jAndrewsQ?! IHadzibabic04j . but such experiments do not require a 
description in terms of atom-number-violating coherent states, as will be discussed 
below. 

As an aside, there is a difference here with experiments in which a single con- 
densate is coherently "cut" into two parts and then allowed to re-interfere. Such 
interference is obtainable without the use of coherent states. As an optical analogue, 
sending a photon Fock state — or even a thermal state — through a Mach-Zehnder 
interferometer demonstrates perfect interference! Experiments in this category in- 
clude |BlochOO] . where interference was demonstrated between matter waves origi- 
nating from spatially separated parts of one trapped condensate, and i Hall9 B| where 
the process of splitting and recombining the condensate was done using two-photon 
pulses, with the two condensates populating degenerate spin states. Only the case 
that the EEC's are prepared separately is of interest in this chapter. 

It is simple to imagine an experiment involving two EEC's that closely follows 
the optical scenario described in Chapter El for photons in cavities. For instance, 
two condensates trapped in separate potential wells may be allowed to slowly tunnel 
through a barrier. Atoms originating from different wells can be rendered indistin- 
guishable by mixing at an appropriate beam sphtter. While a standard description of 
the experiment would utilize interference between coherent condensate fields |^i|e*^ 
and |'?/'2|e*'^, the discussion of the Chapter |21 can be carried over to conclude that 
such a description is not necessary. In fact it is less desirable — it violates atom 
number conservation and invokes the use of the independent and unknowable phases 

^the absolute phase for a single mode is, roughly speaking, the phase difference between the 
amplitudes for some number n and n + 1 particles. 
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9 and 0, which vary from run to run of the experiment and should therefore be in- 
corporated in a quantum mechanical framework by the use of mixed states (leading 
to a description as in Sec. 12.3.11 of Chapter |21). In practice the most striking BEC 
interference patterns are not those which involve leaking of single atoms onto a beam 
splitter and detection in one of only two channels, but rather those in which spatial 
diffraction of the initially independent EEC's occurs, and a spatial interference pat- 
tern is measured in the region of overlap. In what follows the discussion is extended 
to this type of experiment. 

4.2 Spatial interference of independently prepared 
Bose-Einstein condensates 

4.2.1 Modelling the atomic detection process and quantify- 
ing the visibility of the interference patterns 

A simplified illustration of the class of experiment considered in this section is given 
in Fig. 14.11 Two condensates are prepared separately and are trapped close together. 
They are then released simultaneously from their traps and, as they fall under grav- 
ity, the condensates expand ballistically and overlap. A probe beam is shone on the 
region of overlap and the atomic density profile is imaged. High-contrast interfer- 
ence patterns have been observed in experiments performed at MIT |Andrews97j 
and at the ENS |Hadzibabic04j . The latter involved the interference of a total of 



30 condensates, which were prepared in a one dimensional optical lattice, and were 
released and imaged together. Various steps were taken in this experiment to ensure 
that the condensates did not have any prior phase correlation: the lattice period 
was long and the potential barrier high so as to eliminate tunneling between sites; in 
one version of the experiment the lattice was ramped up during the cooling process 
so that the condensates were formed independently; and in another the condensates 
were kept before the start for a "holding time" very much longer than the observed 
decorrelation time. 

It is not necessary to use a coherent state description to predict interference 
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potentials ^ 
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Figure 4.1: A simplified illustration of a spatial interference experiment for two initially 
uncorrelated Bose-Einstein condensates. The condensates are released from their traps, 
expand ballistically, and overlap. The atomic density distribution is imaged optically. 
Interference fringes are observed in the region of overlap. 



between independently prepared EEC's, as was first shown in detail by Javanainen 
and Yoo |Javanainen96bt IYoo97j . They showed that interference patterns emerge 
even if the atom number superselection rule is obeyed exactly, and the condensates 
are initially in atomic Fock states with the same number of atoms. The analysis 
which follows is based on the same simple model employed by Javanainen and Yoo, 
but considers the case that the initial state of each condensate is mixed. It is assumed 
that phase diffusion, the shape of the trapping potential, and edge effects can be 
ignored. The condensates correspond to macroscopic occupation of one-particle 
plane-wave states with opposite momenta ±k. They can therefore be described 
by second quantised fields of the form e^'^^^b±k, where the b±k are annihilation 
operators for each mode. The initial states of the condensates are taken to be 
Poissonian, consistent with a superselection rule for atom number. This choice of 
initial state is suggested by a maximum-entropy principle which can be applied here, 
given the situation of continuous exchange of atoms between condensed atoms and 
a surrounding thermal cloud where only the average size of the condensate is easily 
determined |Jaynes03| . Both condensates are assumed to have the same expected 
atom number N. To generate the pattern of spatial interference, the condensates 
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are assumed to merge over a linear array of atom detectors, with the atoms detected 
one at a time. The combined field operator, i^ = e^^^^bk + e~*^^j6_fc, serves as the 
measurement operator for a detection at position Xj. 

The situation here turns out to be dynamically equivalent to the optical problem 
discussed in Chapter El when the cavity modes are initially Poissonianly populated 
and the second cavity undergoes random phase shifts between detections, for exam- 
ple because of a frequency mismatch (see Sec. I2.1.3[ third paragraph). Inspecting 
the atomic measurement operator t^, it is seen that atomic measurements n/k apart 
are equivalent, and further, that a detection in ^ < a;i < | is equivalent to one in 
< Xi < ^ at a;i — ^ with operator = e^^^^hk — e~*^^'i&_fc. For a mixed state 

poc j d'^ad'^p P{a,p)\a){a\k'S)\P){P\-k, 

the probability density for measurement at Xj with this periodic identification is 
proportional to 

oc j dW/3P(a,/3)(|a|2+ 

On this reduced range every xj is equally probable and the problem can be treated by 
assuming "left" and "right" Kraus operators, Kr^r oc hk+e^'^h^k and Ki^r oc hk—e^'^h-k 
with r taking a random value for each measurement. In fact this equivalence between 
the optical and BEC cases holds regardless of the states involved. 

To formulate a measurement process which can be consistently interpreted in 
terms of probabilities it is necessary to look at the normalisation of the measure- 
ment operators ■il){x). For a general field operator the operator (p\x)(p{x) is 
commonly interpreted as the particle density at the position x. The periodicity 
here of ip'^ {x)ip{x) suggests restricting x to a length of vr/A;. It is easily seen that 
j^^'^ {x)tp{x) = (l) (iVfc + iV_fc j where N±k are the number operators for the 
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condensates with momentum ±k. Hence the operators, 



Kix) : p 



(4.1) 



define a complete set of measurement operators (for more on POVM's and quan- 
tum operations refer for example sections 2.2.6 and 8.2 of INielsenOO] ) . To allow a 
consistent probability interpretation account must be taken of the expected total 
atom number for p, giving a final state |l/tr (iVfc + iV_fc)p ^ K{x)pK'^{x), which 
is subnormalised with trace equal to the probability of the measurement outcome. 
The action of the measurement procedure on the vacuum state with no atoms pop- 
ulating either condensate mode is specified as to resolve the ambiguity in this case 
i.e. the probability is for detecting an atom anywhere in this case. Defining 



/(Xi, X2,--- , Xr) =MK{Xr) ■ ■ ■ K{xi)pK\xi) ■ ■ ■ K\ 



(4.2) 



for a sequence of measurements at positions xi, x^, where M is the normalisation 
factor correcting for the changing expected total atom number, it is easy to check 



that Jq^'' f{xi)dxi = 1 when r = 1. For general r, f^^'' f{xi, ■ ■ ■ ,Xr-i,Xr)dXr = 
/(xi, ■ ■ ■ , Xr-i), and hence /(xi, ■ ■ ■ , Xr) is a well-defined probability density. 

Consider now the effect of the measurement operator K{x) on an arbitrary state 
p for the two condensates (but not the vacuum!) — for example the state after a 
series of earlier atomic detections at different positions. Then the probability for 
detection of an atomic at a position x is given by. 



ir/k 



P{x) 



tr 
1 + 



'k 

i2kx 



p /tr 



tr [hphlA + e-'''^tr [h.kphl] /tr 



hlhk 



b-kh-k ) P 



tr {hkph^_^ is a complex number and tr (^b-kpbl^ its conjugate. Defining 



k P 



tr ( bkpbl 



tr 



k P 



V 



exp-*^. 



(4.3) 
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it follows that the probability density takes the simple form, 

P{x) = l + Vcos{2kx-ip), (4.4) 

where ^ is a positive number between and 1, a fact which is guaranteed by 
the fact that the density operator ip^xYi/j^x) is positive. The result Eq. ()4.4|) was 
identified in |Javanainen96b| IYoo97l but only for the case when the initial states of 
the condensates are number states with the same number of atoms. 

The form of Eq. ()4.4|1 suggests as a measure of the degree of localisation of the 
relative atomic phase parameter of two condensates. On the other hand, in Sec. 12.21 
of Chapter |2l a visibility was defined to quantify the localisation of the relative 
phase for two modes as illustrated in Fig. 12. 4| and this definition can be translated 
to the context of atomic optics. An experiment can be imagined in which one of 
the condensates is subject to a controlled phase shift r and coherently combined 
with the other by a Josephson coupling in analogy with an optical 50 : 50 beam 
splitter. The atom number in one well could then be counted to provide an estimate 
of the intensity /(r) defined in Sec. 12.21 In principle the experiment would have to 
be repeated many times with identical preparations of the condensates to provide 
an accurate estimate of /(t), and repeated again for different values of r so as to 
evaluate the visibility according to equation Eq. ()2.12j) . However there is another 
interpretation. It is straightforward to show mathematically that. 



tr 



(4.5) 



Therefore /(r) is proportional to the probability density P (— ^) for detecting a 
particle at x = — ^. Inspecting Eq. ()4.4|) . V is seen to be the same as the visibility. 

It is also worth noting the connection between the visibility and the method used 
by M0lmer for quantifying the degree of localisation of the relative optical phase in 
his studies of the interference of independent light sources |M0lmer97a| |M0lmer97b| . 
The quantity considered by M0lmer is ?7 = A23? |tr [a"'"6p (t)] }, proportional to the 
real part of the trace quantity, where t is time, a and b are annihilators for the optical 
modes, p is the state of the cavities, and A is a normalisation factor scaling t] within 
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±1. M0lmer restricted his attention to the case of initial number states, in which 
case A is the reciprocal of the number of remaining photons, and he normalised 
the state of the two optical modes. M0lmer assumed a frequency difference between 
optical cavities, giving rise to a rotating phase factor in the governing Hamiltonian — 
h {ujb — uja) b'^b — which causes f] to oscillate. When this oscillation is the principle 
evolution, 77 ^ 23? [Zg^Cv+K-'^^M] = 1/cos[v9+ (cufe-cu^)^], and the visibility is 
seen to govern the magnitude of oscillation of rj. 



4.2.2 Analytical treatment for Poissonian initial states 

This section looks at how the localisation of the relative atomic phase of two con- 
densates induced by a series of measurements by operators defined in Eq. ()4.1|1 can 
be treated analytically. It is supposed that the condensates are initially in Poisso- 
nian states with the same mean number N. The problem can be treated in terms 
of "left" and "right" Kraus operators Ki -^ oc bk — e*'^6_fc and K^^r = bk + e*'^6_fc. To 
understand the characteristic localisation it is sufficient to take half the detections 
at r = and the rest at r = 7r/2, the largest difference possible. There is little 
advantage working as in |M0lmer97al |M0lmer97bl ICastinQTj with a probability den- 



sity for the full measurement record involving information about the precise spatial 
distribution of the atomic detections. The commutativity of Ki r and Kr^r allows 
the process to broken down as convenient. It is assumed that there are M measure- 
ments at each of r = and r = tt/2. The numbers of "left" and "right" counts are 
denoted by /i, ri at r = 0, and I2, r2 at r = 7r/2. 

The 2M measurements cause the initial state with average atom number N for 
each condensate 

Pi = I ^\a){a\,0\P){PU 
where a = y/Ne''^ and P = VWe'^'^, to evolve as 

nr ML/J-'-a f{h T{riT^h /V-'it Tv-nt 7>'«2t T^rat 

'r„l/„l'rJ7J ^.f '.f ^'0 ''O'^^ ''0 '■'0 '.f 



rsUslril/i!' 



-k ' 
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where |C/^^(6', 0)|^= ^^^^ \'^'^^ i^Y^l"^^ |sin (^^)|^' and A = — The peaked func- 
tion Ci^r is famihar from Sec. 12.11 of Chapter |2l and the phase shift r causes a 
translation. 

The general features of the localisation are as in the optical analysis in Sec. 12.3.11 
of Chapter 121 However the effect of the phase shift is to ensure with high probability 
that exactly one value Aq for the relative phase is picked out. This is so even 
when M is small. The phenomenon can be understood by careful inspection of the 
measurement record and with the aid of the asymptotic expressions for C;^r(6',0), 
Eq. ()2.4|) and Eq. ()2.5|) of Chapter [21 Consider now the cases of M = 3, 8 and 15, 
looking at the "likely events" — defined as those with probability greater than a 
equal fraction 1/(M + 1)^. The probabilities of these events total 0.9, 0.8 and 0.8 
respectively. In every case a unique value of Aq is picked out. In very many cases — 
all when M = 3 and half when M = 15 — all the detections are of the same sort, all 



Ki (or all Kr) , at T = OT T = f , or both. In other words at least one component of 

2 2 2 

is of the form |CAf,o| {or |Co,a/| ) which has only one peak 



and a larger spread than otherwise. The product 



2 



in turn 



has only one peak and is highly probable. Other probable events are such that one 



|2 



2 



2 

. In short, 

2 



~ 2 
^l2,r2 



to 



peak of ri(^) r strongly overlaps with one peak of 

I 1 2 

the phase shift of 7r/2 makes it impossible for \C[^~i^^{6 , (f))\ and 
strongly reinforce each other at more than one value for the relative atomic phase. 
In the limit of a large number of detections, 

'de \a){a\k® \ae'^°){ae'^°\.k- (4.6) 

When the relative phase is perfectly defined the atomic detections have a probability 
density of cos^(A;xi — Ao/2) (where xi is the proper position). However, beginning 
from initial states with no relative phase correlation, the value for the relative phase 
localises much faster than it takes for the characteristic spatial interference pattern 
to become well established (a point also made in |Horak99] for the case of initial 



number states). The numerical studies of Javanainen and Yoo, for example, simulate 
interference patterns based on 1000 atomic measurements. However, the dependence 
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on the total number of detections / + r in Eq. ()2.4j) and Eq. ()2.5|1 of Chapter El 
demonstrates that the underlying rate of localisation here is similar to that at either 
of the two values of Aq which evolve when the phase r is fixed. 

Beyond the very first few measurements, the scalar function \C[_^^^Ci^ ^.^ p is well 
estimated by a Gaussian with width between \/2/ v^M and 2/v^M. The state of the 
condensates then takes the form, 



PG 



2n r2iT 



JO 



27r 271 



27rP(A) + (contributions for 27r-periodicity in A) 



X 



Ne 



i(f> 



(4.7) 



where A = — 6* and, 



P(A) = ^ g-(A-Ao)V(2.^) 

(T\/27r 



is a normal distribution with mean Aq and variance a^. The Gaussian width 2(T is 
small compared to 27r (the range of A), and 



2M 



<a < 




(4.8) 



where a total of 2M measurements have been performed. The visibility for a state 
of the form Eq. ()4.7|1 is derived in appendix |Bl and. 



V 



(4.9) 



4.2.3 Numerical simulations for Poissonian initial states 

This section discusses how the emergence of a pattern of spatial interference can be 
simulated numerically, extending the method described in |.Tavanainen96b[ IYoo97j 
to general initial states, and looking specifically at the case of initial Poissonian 
states considered in the previous section. 

The approach used by Javanainen, Yoo and Ruostekoski to simulate the inter- 
ference between condensates initially in number states with the same number 
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of atoms makes key use of the simple form of the probabihty density for atomic 
detection in their measurement model — P{x) = 1 + ycos(2A;x — (p) — and is as 
follows. V is initially 0. V and (p are updated after each atomic detection, and 
a simple method is used to generate the subsequent detection according to P{x), 
assuming a random number generator which supplies a number u on [0, 1] according 
to the uniform probability density. It is easy to show in general that given a random 
variable X, with a probability density function g{x') and a cumulative probability 
function G{x) = J^^ g{x')dx', the random variable G{X) is uniformly distributed. 
Hence to generate the next detection it is sufficient to solve the following equation 
in x: 

r V 

G{x) = / l + Vcos{2kx' -ip)dx' = x + — [sin {2kx-ip) + sm{ip)] = u, (4.10) 
Jo 2k 

which is easy to do numerically in practice. To update V and ip after this detection 
it is sufficient to evaluate P{x) at two points by evaluating the full expression for 
P{x) in terms of the positions of all the previous detections and A^. Since the simple 
form of P{x) in Eq. ()4.4|) applies for arbitrary states, the same method of stochastic 
simulation can also be applied to initial states other than number states. 

In what follows the case that the EEC's are initially in Poissonian states is 
considered in detail. After measurements at positions xi, - ■ ■ ,Xr the initial state 

Pi = I I Wl, (4.11) 



^0 



where a = V Ne'^^, P = V Me^^ and, N and M are the initial expected atom numbers 
for each condensate, evolves to 

P' = ^^l^y K{xr)...K{x{)pjK{x,)\..K{xr)\ (4.12) 

Here K{xj) is the measurement operator for a detection at position Xj (refer Eq. (@^) 
and account has been taken of the total atom number at each detection, p' is sub- 
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normalised such that tr (p') is a probabihty density. It can be shown that 



tr (pO 



{N + M) 



(N + M) 




(4.13) 



where Cm is the coefficient of a"^/?*" in the product H^^^ ( 



easily computed. This expression can be used to update the simulation parameters 
V and (fi after each detection. 

An interesting point concerning the measurement operators Eq. ()4.Hl emerges 



— and it would appear that atom number is not conserved in this analysis. In 
fact more extreme examples are possible. If the two atomic modes are initially 
thermal and both with mean atom number ( A^) , then the total atom number after r 
measurements can be shown to be {N){r + 2). This difficulty can be understood by 
reference to a quantum trajectory type approach which does conserve atom number 
(and as is used in |Cirac96j for example). In such an approach there is, in addition 
to jump operators corresponding to atomic detection, an additional non-Hermitian 
evolution with Hamiltonian ifcond = (^^k^k + b^^k^~k^ where T is some decay 

rate, corresponding the effects of continuous observation of the combined condensate 
conditioned on no atoms being detected. Hence the measurement operators Eq. ()4.1|] 
describe both a detection of an atom at a specific position and a partial measurement 
of the total atom number. 

It is necessary to check whether neglecting the evolution due to -ffcond affects 
the basic features of the localisation process. For the case of initial number states 
|A^)|M) this is straightforward to see. After r atomic detections the state of the 
condensates in a basis of number states is a superposition of states with total atom 
number N + M — r. Therefore the evolution due to -ffcond merely affects the nor- 
malisation of the state. However the case when the condensates begin in mixtures 
of states with different total atom number — and therefore different decay rates — 
is less clear. That the evolution due to -ffcond is not important for the case of initial 
Poissonian states can be understood as follows. A pair of coherent states evolves 



in the explicit calculation of p' above — tr 



blbk + b^_kb-k)p'] = (N + M)tr(p') 
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under -ffcond as 



ie 



Me' 



-(l/2)(Af+M)(l- 



and the conditional evolution is seen not to change the ratio \fWjM of the ampli- 
tudes of the coherent states. It is this ratio, rather than the absolute values of iV 
and M, which is the relevant factor in the probability density after r measurements 
Eq. ()4.13j) (and hence also the parameters V and <y9 in Eq. ()4.4j) ) since 



{N^M''~') / (iV + M)" = 1/ I [1 + {MIN)\ [(AT/M) + 1]"" 



The same is true the scalar factor |c 0)|^ describing the localisation in the relative 
phase A = — 6*: 



P7tr(p') = J^^\c{e,<j>)na){a\ 



and, 



ic(e,0)r 



tr 
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j=0 



Cj Ni M^ e-'^^ 



(4.14) 



The results of repeated simulations for Poissonian initial states with the same 
mean atom number for both condensates are shown in Fig. 14.21 and Fig. 14. HI The 
visibility V for each run is always 0.5 after one detection, but varies significantly 
from run-to-run beyond that, and evolves rapidly to 1 in all cases. In Fig. 14.21 each 
point shows V after a given number of detections, averaged over 5000 runs, and 
the bars mark ±1 standard deviation. Fig. 14. HI serves to verify the prediction made 
at the end of Sec. 14.2.21 which specifies that the scalar function \C {6,(f)) P can be 
well described by a single Gaussian with standard deviation in the range Eq. ()4.8p 
beyond the first few detections. From Eq. ()4.9j) the quantity —2D logg V, where D 
is the total number of atomic detections, is expected to vary between 1 and 2 where 
the result applies. The solid lines mark the predicted range of — 2D logg and 
the dotted curves show the results of the numerical simulation. The central dotted 
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curve corresponds to the average value for —2D logg V and falls within the expected 
bounds. The lower and upper dotted curves are computed for the average visibility 
plus and minus one standard deviation. It may be concluded that the analytical 
prediction is substantially consistent with the simulated results, although the real 
rate of localisation is slightly faster than estimated. 

4.2.4 Comments on the information available in principle 
from a pattern of spatial interference 

This section discusses two open questions concerning the information about the 
relative atomic phase available, at least in principle, from a given pattern of spatial 
interference. The first question is whether full knowledge about the relative phase 
is contained in a precise record of the positions of the atomic measurements, or 
whether some of the phase information is "hidden" , and can only be inferred using 
an additional experimental procedure. 

Specifically the situation imagined here involves an experimentalist who is not 
concerned with the physics of the localisation of the relative atomic phase, and who 
supposes that the EEC's have fixed but unknown phases prior to the experiment. 
It is assumed that the experimentalist is able in principle to fully characterise all 
relevant features of the experiment including: interferometric information - such 
as phase shifts due to path length differences; the values of all control parameters 
determining, for example, the coupling of each condensate out of its potential well. 
The experimentalist is supposed to employ a highly efficient and accurate method 
of atomic detection, and to keep a detailed account of how the condensates are 
prepared. 

Adopting the standard approach of Bayesian statistical analysis, the experimen- 
talist would update the probability density P{A) summarising his or her "knowl- 
edge" of the relative phase distribution after an atomic measurement at position Xr 
according to the rule. 





(4.15) 



/o''/(x^|A')PP"°'- (A') dA' 
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The question is then if there are circumstances where there are substantial differ- 
ences between P (A) and the actual quantum relative phase distribution. Consider 
for example the situation considered in Sec. 14.2.21 for which the condensates are 
known initially to be the same size and uncorrelated. The prior probability density 
is flat. As a pattern of spatial interference is built up, a measurement at a position 
Xj updates the Bayesian probability density according to a factor proportional to 
cos^ {kxj — A/2). After a series of measurements at Xi, ■ ■ ■ , Xr the posterior distri- 
bution is proportional to 11^^]^ cos^ [kxj — A/2). However, this is the same as the 
scalar factor \c{6,(f))\'^ describing the localisation of the relative phase parameter in 
a basis of coherent states, when the condensates are initially Poissonian. The fast 
localisation of the relative phase in this example is mirrored by the rapid narrowing 
of the naive Bayesian phase distribution. 

As another example consider the case in which a pattern of spatial interference 
is obtained for a state localised at two values for the relative phase, 

a){a\ ® [lae^'^-'Xae'^"! + |ae-'^«X«e-'^°l] , 

where |a) is a coherent state with parameter a = Ai/iVe*^, prepared in an atom- 
optical version of the experiment described in Sec. I2.H.1I of Chapter |21 If equal 
numbers of atoms were detected at the two "output ports" during the preparation 
procedure then Aq = vr/2 and the visibility is 0. After the first detection of the 
spatial interference experiment at a position xi, ps evolves to 

p's = cos' (fca; - 7r/4) / — |aXa| ® |ae'"/'X«e'"/'| 

Jo 27r 

+ cos' {kx + 7r/4) / — |aXa| ® \ae-'''/^){ae-'^/% 
Jo 27r 

where ±k is the momenta of the condensates as they overlap. If xi is near either 
^ or ll one or other of the components of ps is picked out and the visibility jumps 
to approximately 1, and the fact that the relative phase is perfectly localised at one 
value is not contained in the interference pattern. However it could be argued that 
the experimentalist who assumes a priori symmetry breaking would not be surprised 
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by this, as the naive Bayesian probabihty distribution, informed by the outcome of 
the preparation procedure, is again the same as for the full quantum state of the 
condensates. 

A second question is whether it is feasible to diagnose the characteristic locali- 
sation of the relative atomic phase by looking at the variation in the points making 
up the pattern of spatial interference. One might expect that whenever the rela- 
tive phase localises as the sequence of atomic detections is made, the pattern is less 
well defined after some number of detections compared to the case when the phases 
are fixed a priori. This effect might be expected to be more pronounced when the 
localisation process is slow. Note that the challenge here is different to the one 
usually considered in quantum parameter estimation, where the aim is to efficiently 
characterise some quantum parameter of a given quantum state by sampling a large 
number of identically prepared systems |D'ArianoOOj . 

Yoo, Ruostekoski and Javanainen compared the case of initial number states 
with atoms in each mode and that of initial coherent states in |Yoo97j . They 
looked at the statistics of the differences X2 — xi, X4 — X3, and so on, for a given 
measurement record Xi,X2, and performed a statistical analysis to determine the 
least number of runs that would be required with initial number states, in order that 
at least half would be statistically significant against the null hypothesis of initial 
coherent states (at the 5% level). The number of repetitions was found to increase 
rapidly with A^ — for A^ = 1, a total of 40 repetitions would be needed, for A^ = 5 
the number is 700 — and too fast to form the basis of a feasible experiment. 

Another approach would be to investigate how the variation of estimates of the 
relative phase changes as increasing numbers of atoms are detected. Given a record 
of atomic measurements at positions Xi, ...,Xr the relative phase could be estimated 
in various ways, for example by maximum-likelihood estimation or by averaging 
single-shot estimates. For the situation considered in Sec. 14.2.21 this would mean 
extremising cos^ {kxj — A/2) in the variable A, and averaging the Aj = 2xj/k 
respectively. However the Cramer-Rao lower bound of classical statistics suggests 
that such a method would not be particularly sensitive to rapid localisation in the 
relative phase. Even for a classical probability distribution where a random variable 
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is phase locked with probabihty density function proportional to cos^ {kx — Ao/2) 
(with Aq and k fixed), the variance of any unbiased estimator of Aq can be shown 
to be at least 1/r for a number r of detections. In conclusion, experiments in which 
the relative phase parameter is expected to localise at a slower rate may be better 
candidates for seeking statistical signatures of the localisation of the relative phase. 

4.3 Non-destructive measurement of the relative 
atomic phase by optical methods 

An interesting direction for further work would explore experimentally feasible ap- 
proaches to probing the localisation of the relative phase between two EEC's, pre- 
pared separately, after an interference experiment involving only a small number of 
atomic detections. It is easy to imagine an experiment in which only small fractions 
of two separated condensates are extracted, coherently combined and subjected to 
a process of atomic detection, for which a pattern of interference would be obtained 
involving only a small number of points, and leaving the remainders of the conden- 
sates in a partially localised state. However demonstrating such partial localisation 
would go beyond any experiment reported to date. The discussion below reviews in 
simple terms a body of work concerned with optical detection of the relative atomic 
phase between condensates. This work treats the issue of the localisation of the 
relative phase briefiy or not at all (an exception being |Ruostekoski97aj ) . but it is 
of interest to identify to several techniques which could be applied to the current 
problem. 

A recent experiment |Saba05j has demonstrated the non-destructive measure- 
ment of the relative phase for two spatially separated and independently prepared 
EEC's, as follows. It is straightforward to show that the density operator 'tp'' {p)^{p), 
for a combined field operator ip in momentum space, exhibits interference when 
the component condensates are spatially separated but phase locked. The method 
used in the experiment to probe this interference is based on |Pitaevskii99j . Two 
counter-propagating lasers are used to continuously extract atoms from the con- 
densates which are trapped throughout the procedure. The lasers interact with the 
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condensates by a process of Bragg scattering (which continuously imparts recoil mo- 
mentum to a fraction of the atoms at high energy and momentum transfer), and 
outcouple them. In addition the Bragg beams overlap throughout a region covering 
both condensates, so that the atoms exit in a single stream, with their origin from 
either condensate indistinguishable. The rate of outcoupling depends sinusoidally 
on the relative phase of the condensates. An energy difference between the conden- 
sates causes the relative phase to change, and so a fringe pattern is observed. For 
each atom extracted a photon is transferred between the Bragg beams (conserving 
momentum), and hence the atomic fringe pattern is mirrored in variation of the 
beam intensities. 

A number of theoretical papers propose schemes to measure the relative phase of 
two condensates spectroscopically, fully or partially. Of particular note is 
|Imamoglu97| which considers a setup with condensates trapped in a two well po- 
tential, and coupled via an excited trap state. The ground state condensates are 
assumed sufficiently far apart to allow them to be driven separately by lasers. In 
one version of the scheme a laser acts as a "local oscillator" and illuminates one of 
the condensates. It is predicted that light is scattered from the other condensate 
with a phase difference (compared to the reference laser) proportional to the one 
between the condensates. The optical phase difference could be measured by optical 
homodyning or heterodyning techniques. This proposal is particularly interesting 
in as much as it suggests that it might be possible to map the relative phase for the 
condensates onto an optical state which would be easier to probe experimentally. 

Three further schemes aim to probe the relative phase of two condensates by 
exploiting features of the light which is scattered when the two condensates occupy 
degenerate electronic ground states, are trapped together, and are coupled via an 
(electronic) Raman system. In each case it is supposed that the Raman system 
is driven off-resonantly by phase-locked, copropagating lasers with different polar- 
isations. In jJavanainen96aj simultaneous amplification of one Raman beam and 
attenuation of the other is identified as a signature of a well defined relative atomic 
phase. |R,uostekoski9 Taj also focus's on intensity variation in the scattered beams, 
and looks at the localisation of the relative phase given initial number states by 
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means of numerical simulation. Finally in |R,uostekoski97bj features of the spec- 
trum of the scattered light are found to depend on the value of the phase difference 
between in the condensates. Another scheme |Savage97| is closely related to these 
three and considers the situation in which excited atoms, prepared independently, 
are launched towards two condensates in overlapping ground states into which they 
can spontaneous decay. Features of the emitted light are found to depend on the 
phase difference of the condensates. 

Finally |Corney9'8| considers theoretically the situation in which Josephson type 
oscillations between two spatially separated EEC's in a double well potential are 
measured by exploiting a dispersive interaction between one of the condensates and 
the electromagnetic mode of a cavity. It is assumed that the cavity mode is off- 
resonantly coupled to an atomic transition, and that the cavity is strongly driven 
by an optical field which leaks out of the opposing end mirror on the time scale of 
the condensate oscillations. The optical signal field and the field which leaks out 
of the cavity are compared using a homodyne measurement scheme. The intended 
interaction between the condensate and the cavity mode has a Hamiltonian operator 
proportional to b'^bc'^c (where b and c are the corresponding annihilation operators), 
and so the optical phase is modulated by the varying atom number in the cavity in 
a manner depending on the relative phase. 
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Figure 4.2: Simulated evolution of the visibility V for the case of initial Poissonian 
states with the same mean number of atoms for both condensates. V is averaged over 
5000 runs. The error bars show ± one standard deviation. 
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Figure 4.3: The quantity —2D logg V, where ^ is a visibility, is plotted against the 
number of atomic detections D. The dotted curves are the results of numerical simulations. 
For the centre one, V is the average visibility after D detections. The lower and upper 
dotted curves correspond to the average visibility plus and minus one standard deviatation 
respectively. The solid lines show the range predicted beyond the first few detections by 
the analysis of Sec. 14.2.21 
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5 



Joint Scattering off Delocalised Particles and 
Localising Relative Positions 



A recent article |R,aun8j by Rau, Dunningham and Burnett (RDB) examines the 
localisation in relative position between two "quantum mirrors" in a Mach-Zehnder 
interferometer (a "rubber cavity"), or between two delocalised free particles, induced 
by the scattering of light. Simple models of the scattering processes are studied 
numerically using a stochastic approach. RDB emphasize the role of entanglement in 
the localisation process. Furthermore they claim that the localised relative positions 
have the properties of classical vector displacements when three or more particles 
are involved, a statement which they substantiate in |Dunningham04 | . This chapter 
presents work published in [CableOS] which continues with the same two models of 
scattering and advances the discussion considerably. 

The two models of scattering in |Rau03j are explained in Sec. 15.11 For the 
first example — that of a rubber cavity — the localisation of the relative position 
between two mirrors closely resembles that discussed in Chapter |21 for the relative 
phase between two optical modes. Differently the localisation of the relative position 
is many valued if the scattering photons have one fixed frequency (with the localised 
positions having separation on the order of the wavelength of the light). 

The second model is for scattering off free particles in two dimensions. It is 
assumed that each particle taken separately acts as a perfect point scatterer, de- 



96 



fleeting the incident light in every direction with equal probability, and with any 
internal excitation of the scatterer relaxing rapidly compared to the time for obser- 
vation. Detection of the photons in the far field renders scattering off two particles 
indistinguishable. Differently from the rubber cavity model the momentum kick on 
each photon detection is variable, causing the relative position of the particles to lo- 
calise at one value. RDB assume in their analysis that the particles are delocalised 
in a one-dimensional region, and it is pointed out here that there are difficulties 
extending the arguments to the case of particles delocalised over a two dimensional 
region. 

The technicalities of analysing the localisation processes in a basis of Gaussian 
states are detailed in Sec. 15.21 A basis of Gaussian states has some advantages 
compared to position eigenstates (which were used in |R,aun3| Dunningham04| ) . In 
particular position eigenstates disperse infinitely fast under free evolution, and an 
analysis in terms of Gaussian states might be expected to run in close analogy with 
the studies of localising relative phase in Chapters El El and E] (which employ coherent 
state representations). However the situation here turns out to be more complicated 
since Gaussian states are not eigenstates of the translation operators corresponding 
to momentum transfer. 

The localisation of the relative position under the second, free particle, model of 
scattering is considered in Sec. 15.31 for the case of particles with equal mass, each 
one initially in a thermal state (with momentum given by the classical Maxwell- 
Boltzmann distribution). In |R,aun3j the initial states for each particle are momen- 
tum eigenstates which is not very realistic. The quantum state of a thermal particle 
has a simple form as a diagonal mixture of Gaussian states (over the position of the 
centre of mass), with spread given by the thermal de Broglie wavelength. Results 
are presented for the cases when the incident light is monochromatic and thermal. 
The situation considered is that of an observer viewing a distant light source — the 
incident light is either forward scattered by the particles into the field of view of the 
observer or is deflected, in which case the light source is observed to dim. Under 
these assumptions the localisation of the relative position between the particles is 
found to be partial, in contrast to the sharp localisation reported in |R,aun8j . For the 
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5.1 Scattering in a rubber cavity and off delocalised free particles 



monochromatic case the pattern of locahsation is characterised by Bessel functions 
of the first kind. 

Finally possible future calculations are suggested in Sec. 15.41 In particular it is 
observed that a substantial simplification of the calculations occurs for the important 
example of scattering off a classical ideal gas at room temperature. Entanglement 
does not play a role in this example, and the free dynamics is straightforward to 
compute. 



5.1 Scattering in a rubber cavity and off delo- 
calised free particles 



Figure 5.1: Photons with momentum k pass through a "rubber cavity" — a Mach- 
Zehnder interferometer in which two of the mirrors are mounted on "quantum springs" and 
are initially delocaUsed along an axis. Two photodetectors monitor the output channels. 

The first model of scattering discussed in |Rau03j — that of a "rubber cavity" — 
is illustrated in Fig. 15.11 (Fig. 1 of |R,aun8j ) . The relative position of two mirrors in 
a Mach-Zehnder interferometer is localised by a series of single photons which pass 
through the device and are detected by photodetectors monitoring the two output 
channels. The Kraus operators corresponding to the two possible outcomes for each 
input photon can be derived as follows. Letting pmirrors denote the initial state of 
the two mirrors, labelling the two optical modes a and b, and working in a Fock 
basis for the optical modes so that the state of the input photon is |10)afe, the action 






k 
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5.1 Scattering in a rubber cavity and off delocalised free particles 



of the first beam splitter is given by: 

Pmirrors 

ab 

^ Pcombined Pmirrors ^ 



|io),,+|oi),,V(ioL.+(oi| 



V2 J\ V2 



lab 



The action of the hght in mode a is to give a momentum kick \/2kNa to the de- 
locahsed mirror along its path (where Na is the number operator and k is the 
momentum of the input photon), and is described by the operator ey:.Y>{i\/2kNaS^ 
(where x is the position operator for the mirror). Similarly for mode h the operator 
is exp ^i\/2fc-^fe^y After the light in the interferometer interacts with the mirrors 
the state of the system is, 

exp (^iV2kNby^ exp (^iV2kNaX^ Pcombined exp (^—i\/2kNaX^ exp 

The second beam splitter, identical to the first, transforms the state of the system 
according to. 



exp{iV2kx) + exp{iV2ky) |01) 



ab 



V2 

exp(^i^/2kx) — exp(iv^fcy) 



exp(i\/2kx) + exp(z\/2fcy) 



m)ab- 



The photon is detected at one of the two output channels, and the Kraus operators 
are proportional to, 

exp(^iy/2kx^ ± exp(^iy/2ky^ . (5.1) 

The localisation of the relative position in this example resembles that of relative 
optical phase in Chapter |2l In fact the (numerically produced) Fig. 2 of |R,au03j 
is essentially identical to the (analytic) Fig. 12.21 in Sec. 12.1.11 of Chapter |21 The 
action of the Kraus operators in a basis of position eigenstates is analogous to that 
of the optical operators Ki and Kr on optical coherent states, although in the latter 
case only the basis is overcomplete. Differently from the optical case the pattern of 
relative spatial localisation which emerges has a periodicity of 7i\^/k and extends 
throughout the region where the mirrors are at the start. 
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Figure 5.2: Plane wave photons with momentum k scatter off two free particles, delo- 
calised in a region of length R, and are either deflected at an angle 6 or continue in the 
forward direction. The observer can "see" photons which forward scatter or which are 
deflected only by a small amount. 

The second example considered by RDB is the localisation in relative position of 
two free particles boxed in a one dimensional region due to the scattering of plane 
wave photons. Fig. 15.21 illustrates the situation. The Kraus operators are derived 
as follows. It is supposed that each particle taken by itself acts as a perfect point 
scatterer, scattering as S-waves with certainty. The scattered photons are detected 
in the far field at some angle 9 of deflection. This simple scattering cannot yield 
information about the position of the particle (see |Heisenberg49| ). Rather each 
event imparts a variable momentum kick A; sin 6' with operator exp(ifc sin where 
k is the momentum of the incident photon, assumed to approach perpendicularly. 
For two particles, and making the assumption that the scattering off each particle 
is indistinguishable, the Kraus operators take the form 



where is some phase factor to be determined. These operators are seen to form a 
density. The condition that the two particles should act as a single point scatterer 
when together sets = 0. In addition, unitarity implies the possibility for forward 
scattering and the corresponding Kraus operator is 



(there is just one). These forward scattering events do contribute to the localising 
process and cannot be ignored. 




(5.2) 




(5.3) 
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5.2 Action of the scattering processes in a basis of Gaussian states 

It could be asked why the two particles are not considered to be delocalised in 
a two dimensional region. In fact there is a complication in this case. Following 
the same arguments as above it is easy to write down the relevant Kraus operators. 
Resolving vectors in Cartesian components with the "x-direction" aligned along the 
direction of propagation of the incident plane wave: 



Ks = Y 2^ 6xp (iAk ■ m) cos ■ f ) 

where the vector position operator for the j^^ particle is (x^, y^), f = {x2 — Xi, y2 — yi), 
^ _ ^ xi+x2 ^ yi+y2 "j g^j^^ ^Y^e momentum kick is Ak = k{l — cos 6', — sinO). The prob- 
lem is this — the operators depend on the vector f rather than just its magnitude 
|f|. The localisation process potentially resolves an absolute orientation as well 
as a relative position, with the fixed direction of the incident photons providing a 
reference. 

The key difference between the rubber-cavity and free particle scattering models, 
as demonstrated in the numerical studies of |Rau03j . is that the changeable momen- 



tum kick of the latter localises a single value for the relative position rather than a 
periodic array. Localisation at any value within the initial region is possible. Notice 
however that the symmetry about remains — the scattering process localises the 
absolute value of the relative position, giving two values for the relative displace- 
ment whenever the initial conditions allow for it. Comparison should be made with 
Chapter E] where translations rather than frequency shifts prevent multiple values for 
the localisation of the relative parameter, and the multiplicity is eliminated entirely. 

5.2 Action of the scattering processes in a basis 
of Gaussian states 

In what follows the models introduced by RDB are investigated further. For this a 
basis of Gaussian states has some advantages over position eigenstates — position 



101 



5.2 Action of the scattering processes in a basis of Gaussian states 

eigenstates are dynamically fragile in as much as they disperse infinitely rapidly 
under free evolution, and using Gaussian states facilitates analogy with the locali- 
sation of relative phase discussed in earlier chapters. However the case of localising 
relative position turns out to be more technically complex since Gaussian states are 
not eigenstates of translation operators. This section takes a diversion to set out 
the technicalities of working with Gaussian states. 

The notation \ipk,a,d) is adopted for a Gaussian state with mean momentum k, 
mean position a and spatial spread parameter d, 

/oo 
dxe^H^-'^^GaAx) k) , (5.4) 
-oo 

where Ga,d denotes a Gaussian probability distribution with mean a and spread 



Consider now the effect of a sequence of localising Kraus operators, K{xi,X2) — 
Kij^{xi,X2)--.Ki-^{xi,X2), acting on an arbitrary basis state for the two particles. 
Rather than explicit localisation in the relative mean position parameter 02 — Oi, 
each basis state evolves to a superposition as follows. 



\'tpki,ai,d)®\i^k2,a2,d) (^) J dXidX2dpidp2 

X K{xi,X2)e'P'^"'^e'P'^'''-''^^ 

X |V'pi+fci,ai,d)®|V'p2+fc2,a2,d) (5.5) 

where K{xi,X2) is evaluated at number values and, unless otherwise specified, the 

integral is /^•••/^. This can be understood more simply. Writing K{xi,X2) = 

g+<' ^2 ^'Ce( ^^2^^ ), a product of translation to the centre of mass (where | is the 

cumulative momentum kick) and the localising function Cq{ ^^~^^ ), and assuming 

that the Fourier transform Ce{p) — J dze~^P^CQ{z) may be defined, the final state 

^To simplify the notation the standard deviation d is suppressed when it is constant in consec- 
utive steps. 
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5.2 Action of the scattering processes in a basis of Gaussian states 



of the particles is proportional to, 



|-§+fcl,ai 



| + f+A;2,a2 



(5.6) 



It is seen that the basis states evolve to an increasingly flat superposition over relative 
mean momentum with the phase terms defining the location of the relative spatial 
maxima. In particular Cq{p) ~ e~*2^° when K[xi,X2) enforces sharp localisation 
of the relative position at a single value Aq. 

The picture can be further clarified by changing basis, regarding the Hilbert 
space as a tensor product of spaces for the centre of mass and the relative position 
rather than spaces for each particle. 



Rel' 



The basis states considered above can be rewritten as another product of Gaussian 
states. 



|V'fei,ai,<i)®|'0fc2,a2,<i) OC 



ai+a.2 d 



2 '%/2/COM 



fc2 — fci 



2 ' V2 /Rel 



In the new notation the final state, after K{xi,X2) has acted enforcing sharp local- 
isation at Aq, is proportional to. 



/Cl+fc2+(„ 2 'V2/COM J 



^jf (a2-ai-Ao) 



(5.7) 



The centre of mass component is merely translated. The localisation in the relative 
component is best seen by comparison with the following identity, expanding an 
arbitrary position eigenstate in terms of Gaussians: 

j dpe'^'^'^-^) \i^k+p,a,d) OC e'^^''-^^Ga,d{X) \X) . 



If X is far from a the norm vanishes. 
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5.3 Localising thermal particles with monochro- 
matic and thermal light 

This section looks at the locahsation of relative position for two particles as might 
occur in nature. Rather than the pure momentum states chosen by RDB for initial 
states, the localisation between two thermal particles is considered here. It is sup- 
posed that the two particles have the same mass m and temperature T. At first it 
will be assumed that the localising process can pick out one value Aq for the relative 
position. 

A thermal state for a single particle is given by a mixture of momentum eigen- 
states weighted according to the classical Maxwell-Boltzmann distribution, and can 
be expressed in terms of Gaussian states in a simple diagonal form, 

/^^^y^^exp(-^)|p)(p| = l^l^PoMoA- (5-8) 



The spatial spread parameter of the Gaussian states is given by = y 2mkBT ("^here 
ks is Boltzmann's constant) and may be identified as the thermal de Broglie wave- 
length (the average de Broglie wavelength for an ideal gas at temperature T). When 
the particle is heavy and hot the Gaussian states approximate position eigenstates. 
To attain normalisable states the infinite limits in Eq. ()5.8|) are dropped and the 
initial states are assumed to be delocalised over a finite region R. For two particles 
together. 



pl(X j I daida2 |^0,ai)(V'0,ai I ® \^0,a2){^0,a2\ ■ (5-9) 
JR JR 

Under the action of the localising process K{xi,X2) the initial state pj is trans- 
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5.3 Localising thermal particles with monochromatic and thermal light 
formed as follows, staying with the notation introduced in Sec. 15 .2^ 

Jr Jr J 

X ^-i{Pii^l-o.l)+P2{x2-a2)}^i{p[(x[-a[)+p'2ixi2-a'2)} 

X K{xi,X2)K{x[,X2y 

X \i^p,,a^,d){'^p[,ai,d\ ® IV^P2,a2,d)(V^p!,,a2,d| , (5-10) 

where d'^x^^ abbreviates dxidx[ etc. When the K{xi,X2) operators enforce sharp 
localisation — at Aq say — the final state takes the simple form, 

Pf oc [ [ [ rfairfa2c/V^'^et(f-P')('^2-'^i-^o) 

JrJrJ 

X ^|„|,,,)(^|_^,,J® (5.11) 

The centre of mass and relative positions remain unentangled throughout the lo- 
calising process as is clear for example from Eq. ()5.7p . and as would certainly be 
expected. The two particles however evolve from being separable to being highly 
entangled. This contrasts to the localisation of relative optical phase for two ini- 
tially Poissonian or thermal states, which remain separable despite the emergence 
of strong correlation between them (as discussed in Sec. l2.3l of Chapter 12)). 

The localisation in relative position induced by the scattering operators Eq. ()5.2j) 
and Eq. ()5.H|1 . describing the general case of light scattering off two free particles, will 
now be treated in detail. Differently from jRau03j it is supposed here that there is no 
access to a detailed record of every scattering event. Rather it is assumed that there 
are two types of measurement outcome: "forward scattering" where the incident 
photon continues without scattering or is scattered into a small angle between — e 
and e; and "deflection" where the photon is scattered outside of this range and the 
light source dims. It is necessary then to mix over all possible events constituting 
each of these measurement outcomes. This is a reasonable model for a real observer 
monitoring light from a distant source scattering off two particles, who only has a 
limited field of view and cannot measure the angle of deflection. 
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5.3 Localising thermal particles with monochromatic and thermal light 

As previously the initial states for the particles are supposed to be thermal and 
to be delocalised in a region R. Changing basis to separate out the centre of mass 
and relative components, the initial state of the particles is, 



Pi oc 



i+a2 



V^n °l+''2 d )( Ipr. ai + g2 d 
2 'v^/\ ^' 2 'v^ 



COM 



PReh 



where. 



PRelCX 



^0, 



2 ' V2 



a2-ai d 
^' 2 ■ 



Rel 



Having changed integration variables from oi and 02 to si±22 g^,^^ 22_2J.^ i^iower and 
denote the lower and upper limits of the inner integral which corresponds to 



-^upper 

the relative component of the two particle state. Liower and Lupper depend on the 
outer integration variable , which in effect ensures that the particles remain 
within the original region R. 

Consider first the case that the incident light takes the form of single photons 
at one frequency. With S deflection and F forward-scattering events p^ei evolves to 
p'j^gi as follows, tracing out the centre of mass component at the end,^ 



PRel ■ 




' °l+a2 



' ai+Q2 ' 



(X 2 2 2 ^2~^l d \ 2 y 

2 ' \/2 



G a 



2-°l d \ 
2 ' V2 



2tt 



c/^sin^ + J de cos {!^^) cos 

s 



r^r^c fks'm8r\ fks'm9r'\ \r\/r' 

df^cos j COS j I2X2" 



I Rel" 



(5.12) 



A typical pattern of localisation is shown in Fig. 15.31 which plots the probability 
density P{y — x) oc (^^IPrciI^^) ^'^^ different ratios of "deflection" and "forward- 
scattering" events, where x and y are the (precise) positions of each particle; prior 
to the scattering process P{y — x) is uniformly distributed. In contrast to the sharp 
localisation reported in |Rau03j . limited knowledge of the scattering record means 
that the localisation of the relative position is only partial even after many photons 
have been scattered. The localisation takes the form of complex interference patterns 

^The momentum kick imparted to the centre of mass depends on the angle of scattering but 
the hnearity of the partial trace procedure ensures that this causes no additional complication. 
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Figure 5.3: Probability densities between Liowcr -^uppcr for the relative separation of 
two free thermal particles after 5 photons, each with momentum k = 5, have scattered off 
them, either being deflected into some large angle or continuing in the forward direction. 

The spatial spread parameter d = ■\J 2mkBT '^^ (units 27r/A;). 



rather than sharp peaks. 

As would be expected the degree of localisation for different outcomes is found 
to be insensitive to the precise value of the small parameter e describing the narrow 
range of angles visible to the observer (as are the associated probabilities). Taking 
e = the interference patterns are characterized by Bessel functions of the first kind 
PRc\{y — a;) ~ [1 — Jo{k{y — x)yf [1 + Jo{k{y — x)yf . The localisation is symmetric 
about the origin. Sharp localisation at one specific value — y — x — is possible 
with small probability and occurs when every photon is deflected. 

Consider flnally the case that the particles are initially in thermal states and that 
the scattering fight is also thermal. Each incident wave packet is then described by 
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Figure 5.4: Probability densities for the relative separation of two free thermal particles 
after 5 thermal wave packets have scattered off them, either being deflected into some 
large angle or continuing in the forward direction. The momentum parameter k = 5 and 
the spatial spread parameter d = 0.2 (units 27r/A;). The thermal wave packets have mean 
photon number n = 5. 



the mixture p = X^n [^"/(-'- + ^)"''"^] I'^X^I (where \n) denotes an n photon Fock 
state). The scattering operators Eq. ()5.2|) and Eq. ()5.3|) must be modified here — 
the fixed momentum kick k sin 6 is replaced by the operator Nk sin 6 where N is the 
number operator for the optical mode. Scattering of a single thermal wave packet 
leads to a variable photon count, all detected at a single angle 9 of deflection. A 
typical pattern of localisation is shown in Fig. l5.4l which plots the probability density 
P{y — x) after five thermal wavepackets have scattered; the results are not sensitive 
to the precise value of the small parameter e. 
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5.4 Possible extensions 

This section discusses some possible extensions of the calculations in this chapter. 
The scattering of light off the particles of an ideal gas at everyday temperatures is of 
particular interest, and it turns out that a substantial simplification of the previous 
calculations is possible for this situation. Consider first the definition Eq. ()5.4|1 of 
Sec. 15.21 where \ipk,a,d) denotes a minimum uncertainty Gaussian basis state with 
expected momentum k, expected position a and spread d. The action of the mo- 
mentum kick operator exp{ipx) (transferring momentum p) on these basis states 
has a simple form whenever the length scales 1/p and 1/k are large compared to d, 
exp {ipx) \il'k,a,d) €,'^^°'\ipk,a,d) ■ Now for a particle of a classical Maxwell-Boltzmann 
gas, the state of the particle is mixed and may be represented diagonally in terms of 
Gaussian basis states as in Eq. ()5.8|1 . with zero mean momentum and a spread given 
by the thermal de Broglie wavelength. The thermal de Broglie wavelength is very 
much smaller than the (classically derived) average distance between particles and 
the wavelength of light at optical frequencies (the thermal de Broglie wavelength is 
for example 0.02nm for molecular nitrogen at 300i^). Hence the simple form of the 
action of the momentum kick operator above would be expected to apply for the 
particles of such an ideal gas for parameter values relevant to the scattering of light 
from an everyday light source. 

Two initially uncorrelated particles of an ideal gas would be expected to evolve 
under the action of the scattering processes given by Eq. (jS.ll) . or Eq. ()5.2|) and 
Eq. ()5.3|) . to good approximation as follows, 

daida2 |V^O,ai)(V^O,ai | ® |^/'0,a2)(^0,a2 I 

daida2\K {ai,a2) \^ |^/'o,ai)(^o,ai | ® |V^o,a2)(^o,a2l , (5.13) 

(continuing with the notation of Sec. 15. 3j) . Here the localisation of the relative 
position is contained in the scalar function |fr (01,02) P, and the analysis has close 
analogy with that in Sec. 12. 3.11 of Chapter [21 concerning the localisation of relative 
optical phase for the case of initial Poissonian states. A full analysis would need to 
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check carefully the regimes in which the result Eq. ()5.1H|1 holds. 

Two simple features of the localisation in Eq. ()5.13|) should be noted. First the 
particles stay unentangled throughout the localising process. However in jRau03| 
|Dunningham04| RDB look at examples with pure initial states for which the par- 
ticles do evolve to highly entangled states, and they suggest that entanglement is 
central to scattering-induced localisation of relative positions in general. This is 
seen here not to be true in an important example. Second the free evolution un- 
der the Schrodinger equation of the localised particles has a simple form, shedding 
some light on the robustness of the localised states. For the mixed diagonal state 
Eq. ()5.13p it is easy to see that the free unitary evolution is simply given by replacing 
the Gaussian basis states with the states to which they would evolve separately i.e. 
\4'k,a,d) — ^ exp (^—iHot^ iV^fc.a.d), whcrc Hq is the kinetic energy operator. The free 
evolution of the one dimensional Gaussian state is well known Gohen-Tannoudji77| . 
A Gaussian state of the form \ipk,a,d) remains Gaussian at all times with fixed mo- 
mentum dispersion and steadily increasing spatial dispersion. The evolved state is 
not a minimum uncertainty state and is not of the specific form of the Gaussian 
basis states as defined in Eq. ()5.4|) . Two particles, initially in localised states of the 



form Eq. ()5.13|) . will each have a spatial dispersion of 2^/1/ {2mkBT)^Jl + Ak'^T'^t'^ 
at a time t > after the start of the free evolution (where m is the particle mass, T 
is the temperature and ks is the Boltzmann constant). 

In another direction, the model of free particle scattering treated by RDB and 
in this chapter is rather simple — it is assumed that for a single scatterer incident 
photons are scattered with certainty and with equal probability in every direction — 
and it is of interest to study more realistic models relevant to the scattering of light 
from an everyday source off air or dust. For example, two types of elastic scattering 
dominate the scattering of sunlight off the atmosphere — Rayleigh scattering off 
molecules and Mie scattering off particles larger than a wavelength of the light. For 
Rayleigh scattering the scattered intensity is proportional to (1 -|- cos^ 0)/^^, where 
6 is the angle with the forward direction and is A the wavelength of the light. The 
dependence of the total cross section on A might have implications for the scattering 
of thermal light which has components with different energies. Mie scattering is not 
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5.4 Possible extensions 



strongly dependent on the wavelength of light, but the scattered intensity has an 
irregular distribution with strong scattering in the forward direction. Interestingly 
the patterns of localisation discussed in |Rau03j and in this chapter are characterised 
by the wavelength of the light but for Mie scattering the scatterers are already of 
this order. 



Ill 



6 

Outlook 



This final chapter suggests several ways in which the ideas and methods developed in 
this thesis could be applied to new problems. It is clear that the "modus operandi" 
for Chapters 121 through toElcan be extended, and also applied to many more physical 
systems, and this is discussed first. Many of the wider issues arising in this thesis 
in the contexts of quantum optical processes and Bose-Einstein condensates, are 
also relevant to circuits with Bardeen-Cooper-Schrieffer superconductor components. 
The coherent coupling of bulk superconductors in close proximity by a mechanically 
oscillating superconducting grain has been discussed theoretically, and provides a 
concrete system where localisation of the relative superconducting order parameter 
could be investigated. There are a number of "big debates" with connections to the 
thesis topic. One of these concerns different mathematical descriptions of quantum 
phase measurements, and compares the features of the various approaches. Another 
is why and how time should be treated as a quantum variable. Here processes of 
localisation are relevant to synchronising "quantum clocks". Finally, there is the 
question of the extent to which the theory of decoherence fully accounts for the 
emergence of classicality in quantum systems. States with a well-defined relative 
correlation, of the type whose preparation is discussed in this thesis, have potential 
application here as pointer states — states which are comparatively long lived with 
respect to coupling to an environment. 
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A "modus operandi" for analysing processes which localise some 
relationally defined degree of freedom in different physical systems 

It is clear that the methodology set out in Chapters El through to El can be easily 
extended. The key ideas and methods discussed in this thesis are also relevant to 
many more physical systems. Consider, for example, the localisation of a relative 
angle degree of freedom defined by two spin systems. This problem is relevant 
to proposals and experiments concerned with using measurements on spin systems 
as the elemental operations in a quantum computer. Experimental detection of 
the magnetic force between a ferromagnetic "tip" and a single electron spin has 
recently been reported, using the technique of magnetic resonance force microscopy 
|Rugar04l . For this system the ferromagnet might be treated theoretically as a 
"quantum gyroscope", serving as a reference for the spin-1/2 particle. Angular 
momentum coherent states have several properties analagous to those of optical 
Glauber coherent states (as explained in jPeres95j for example), and might provide 
a mathematically convenient representation along the lines of the analysis in this 
thesis. 

It should be pointed out that there are also many systems where it is of interest 
to study some dynamical process localising a relationally defined parameter, which 
are quite unlike the measurement-based processes studied in this thesis. To pick just 
one example consider an electron and a proton, initially well separated and "free", 
evolving under the Coulomb potential to a stationary state of a hydrogen atom with 
a well-defined particle separation. Taking a simple model the stationary states here 
could just be taken to be those of the time-independent Schrodinger equation for a 
"simple" hydrogen atom (as detailed in introductory texts). 

Localising relative superconductor phase 

The Bose condensation of Cooper pairs (weakly bound pairs of electrons) to a 
macroscopically occupied ground state plays a key role in the phenomenon of su- 
perconductivity. A superconductor is then typically assigned an order parameter 
corresponding to a definite, absolute, phase. However, a line of thinking similar to 
that in Chapter IH concerning atomic Bose-Einstein condensates asks the question 
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of whether this description has possible shortcomings, in particular with respect to 
careful considerations of charge conservation. The coherent coupling of two bulk su- 
perconductors in close proximity by a superconducting grain shuttling mechanically 
between them has been discussed theoretically IGorelikOTI IIsacssonn2] , and provides 
a concrete example where localisation of the relative superconducting phase could be 
investigated in detail. In this system a moveable nanometer sized superconducting 
grain acts as a Cooper pair box. The grain behaves like a two level system defined 
by states different by twice the fundamental electron charge (i.e. by occupation by 
a single Cooper pair). It is able to interact with the bulk superconductors when 
touching them via a suitable gate mechanism. In addition, the bulk superconduc- 
tors can be connected electronically by a Josephson junction. The magnitude of 
the current through the junction depends sinusoidally on the value for the relative 
phase. Future investigations might aim to clarify the roles of the mechanical and 
electronic couplings in establishing well-defined relative phases between the bulk su- 
perconductors, identify the most appropriate description of the initial states of the 
superconductors (including the possibility of mixed states), and explore analogies 
with localising optical and atomic relative phases. 

Quantum phase operators and distribution functions 

A variety of methods exist for treating phase variables in quantum mechanics, 
such as Pegg-Barnett phase operators and Wigner functions. There has been much 
discussion of the features of these different approaches (see for example |Buzek92] 
Garraway92| ). The POVM (positive operator- valued measure) operators a±b (where 
a and b are annihilation operators for two simple harmonic oscillator modes), or 
a + (where the phase shift ^ takes values on a continuous range), summarise the 
measurement processes in Chapters El through to|3J and constitute natural measure- 
ment operators for the relative phase between two modes. It would be valuable to 
study how the "relative number" variable defined by the oscillator modes changes as 
the corresponding relative phase becomes well-defined under the action of these mea- 
surement operators. The situation becomes especially interesting when the initial 
state has a larger intensity in one mode. 
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Synchronising quantum clocks 

Another topic which is related to this thesis is that of "quantum clocks" . Time 
in quantum mechanics is commonly and implicitly taken to be a classical variable 
— "coordinate time". This is the case in the Schrodinger equation for example. 
However, an extensive literature exists which attempts to treat time in terms of 
measurements on systems of quantum resources — spin systems, or systems of sim- 
ple harmonic oscillators for example. The motivation here is to deal with various 
conceptual difficulties associated with using coordinate time in quantum mechanics, 
such are manifested in the well-known Wheeler-DeWitt equation. If time is treated 
quantumly, dynamics that would traditionally be described as a unitary process 
must now be treated probabilistically (as is done in the conditional probability in- 
terpretation of Page and Wootters |Page83 for example). The issue of the dynamical 



localisation of relative degrees of freedom is relevant here for problems concerning 
the synchronisation of quantum clocks. In fact the localisation of relative phase 
between optical modes initially in Poissonian states, as analysed in Chapters |21 and 
ini is immediately applicable — for a laser constitutes a simple, concrete, example 
of a quantum clock (a matter treated at length in |Wiseman04j ) . 



Pointer states for decohering systems 

Another possible direction for this programme of research is an investigation of 
the relational aspects of the theory of decoherence. The situation here is of a focus 
system interacting with a reservoir (an environment). The two are assumed not 
to be correlated at the start. After some time the focus system and the reservoir 
evolve to become highly correlated. In general, it is not possible to keep track of 
the exact state of an environment, and the state of reservoir should be traced out. 
This yields the focus system in a highly mixed state. Overall, the key process of 
decoherence is the damping of the off-diagonal coherences of the focus system in a 
special basis of "einselected pointer states" . An initial state of the form, p^^'^^f = 
'^ijCiC* \4>i){(pj \ , where the are pointer states, evolves to the diagonal form 
pfocus'^^'^'^ — Einselected pointer states have special properties — 

they are particularly stable against interaction with the environment, and they lose 
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their purity only slowly. However, for a given system it is difficult in general to 
determine a basis of pointer states. The existence and completeness of such states 
is not guaranteed, and they can be hard to identify. 

States with a well localised value of a relative parameter, analagous to those 
whose preparation is explained in the thesis, are proposed here as natural can- 
didates for pointer states in a more relational approach to decoherence theory. As 
demonstrated in this thesis, sharp and rapid localisation of a relationally defined de- 
gree of freedom is possible in some systems even when the initial states are mixed, 
and happens as readily for an apparatus with small instabilities or asymmetries 
as in the ideal case. This suggests that such dynamical localisation is not a frag- 
ile phenomenon that can only be observed in highly controlled experiments, but is 
something of relevance to processes occurring in nature. Examples of processes of 
localisation acting on highly asymmetric initial states are particularly relevant here. 

However even after adopting a more relational approach to the study of decoher- 
ence processes, the problem often raised of why definite events are actually observed 
in the real world (see for example the debate in |Zurek9l] [RepliesOSl ) is likely to 
remain. In particular, the correct interpretation of the final mixed state pfocns'^'^'^'^ ^ 
diagonal in a basis of pointer states, remains to be addressed, pfocus'^^'^'^ is an im- 
proper mixture and cannot be interpreted as describing a situation where one of the 
outcomes "actualises" according to the probability distribution It remains an 
open question as to whether or not careful consideration of relationism in quantum 
mechanics has anything to say on this issue. 
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A 

Derivation of the visibilities for Poissonian 
and thermal initial states 



A.l Derivation of the visibility for Poissonian ini- 
tial states 

The initial state of the cavity fields, a product of two Poissonian states both with 
average photon number N, 



I 



dOdcf) I , , I 
a) {a\ (g) 



47r2 

where a — y/Ne^^ and /3 = VNc'^'^, acquires a factor 

^ i\ I v^« + ye/3\ /, -v/ea + v/e/9\ ^. 
0) = -^^) [l (A.1) 

extracting the I and r photon components of the coherent states under the canonical 
localising process, in which a fraction eN leaks out of each cavity and, / and r 
photons are detected at the left and right detectors respectively. \Cir\ is a function 



of A = — 6* and is peaked at ±Ao given by 2 arccos ^/r/r + /. The final state is 
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A.l Derivation of the visibility for Poissonian initial states 



then, 



^r+l 



-2eN 



a + P 



V2 



2r 



-a + p 



V2 



21 



where a' — \/l — ea and (5' — \/l — e(5. 



A. 1.1 Calculation of the probabilities of different measure- 
ment records 

A probabihty can be calculated using, 



/ 



^a^j. 2r^-^ ■ 2i(t>-0 , r(r + 0.5)r(/ + 0.5) 
^ 2 2 r(r + i + l) 



Pi,r{e,N) = Hp' 



(2eiV)'-+'^_2,jvr(r + 0.5)r(/ + 0.5) 



r\l\ 



7rr(r + / + 1) 



A. 1.2 Calculation of the visibility after a sequence of detec- 
tions 

The visibility of p' is computed as follows. The second mode undergoes a variable 
phase shift of r and both modes are then combined at a 50 : 50 beamsplitter, p' 
goes to p" according to 

\a') {a'\ I/?') {I3'\ ^ 



^/2 



And an intensity at the left detector is then given by 



7(t) = tr(a^ap") 



oc 



dOdct) 



a + P 



2r 



'a + P 
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a' + P'e'^ 
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A.J Derivation of the visibility for Poissonian initial states 

where the constant of proportionahty is of no interest (it divides out when computing 
the visibihty). Expanding the last term of the integrand, 



a' + 



72 



2(1 -e)A^cos 



2,A + r, 



where A = (p — 6. The expression for J(r) may be simphfied. 

Kr) oc 



cos-(y)sin^'(y) 

COS ( — ) COS - 

,A, ,A, r r 
—2 cos( — ) sin( — ) cos — sin — 
^2^ ^2^ 2 2 

+ sm (— )sm - 



The first and last contributions can be resolved in terms of Gamma functions, as 
for the probability above, and the second term evaluates to 0. So, 

I(t) oc r cos^ — h / sin^ — I — 
^ ^ 2 2 2 

and extremising at r = and r = tt, 

^ (-^max -^min) / (-^max ~l~ -^min) 

\r — l\ 



r + l + 1 



A. 1.3 Revised calculation of the visibilities for constituent 
components localised at one value of the relative phase 

The localising scalar Ci^r{0, 4') is the same as for the case as initial Fock states with 
the same number, illustrated in Fig. 1221 It is clear that p' may be considered a sum 
of two separate components, one localised at +Ao with the relative phase parameter 
A/2 = varying on to 7r/2, and another at — Aq with A/2 restricted to — 7r/2 
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A.J Derivation of the visibility for Poissonian initial states 



to 0: 



+ 



TirW. 

^r+l 



-2eN 



-2tN 



d 



a + P 



V2 

a + (3 



V2 
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V2 

-a + (3 
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\a'){a'm(3'){(3'' 



V2 
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|a>'|®|/3')(/?" 



In what follows the visibility is computed for the +Ao component separately. The 
computation proceeds as in Sec. IA.1.21 Differently from this previous calculation, 
the reduced range of the integration variable A/2 means that the intensity /(r) has 
an additional contribution (the "cross term" no longer drops out). In detail, 



I(t) oc cos ( — 



tt/2 



dx cos^*"^^ (x) sin^' (x) 







- 2 cos (^) sing 

7r/2 



7r/2 



dx cos^'-+^ (x) sin^'+i 



+ smM2 



r 



dx cos^*" (x) sin^^'"*"^^ (x) 



where the integrands are all positive on the to 7r/2 quadrant. Evaluating the 
integrals in terms of Gamma functions: 

i(t) oc cos" (-)r(r + ^ir(/ + i 



2 COS (^) sin (I)r(r + l)r(/ + l) 



+ sin^ - r r + - r / + - 



Computing the maximum intensity, /max, by setting r = — Aq, and the minimum 
intensity, Imm, by setting r = — Aq — tt, the visibility can be computed: 



^ (/max /min) / (/max ~l~ /n 



(r - if + |4v^r (r + 1) r (/ + 1) / [r (r + 1) r (/ + 1)] } 



(r + /) (r + / + 1) 
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A. 2 Derivation of the visibility for thermal initial states 



A. 2 Derivation of the visibility for thermal initial 
states 



These calculations follow a similar line to the Poissonian case above. The initial 
state of the two cavity fields, a product of two thermal states with the same average 
photon number N, 
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Nn 



d^acfp exp— 



\al+\Pl 
N 
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where a = ^/ne and (5=\/me^'''^ acquires a factor 



V V2 /\ 71 /' 



extracting the / and r photon components of the coherent states under the canonical 
locahsing process. 
The final state is 
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d^ad'^P exp 
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where a' = Vl — and P' — y/1 — ep. This expression may be simplified using the 
parallelogram rule, 
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A. 2 Derivation of the visibility for thermal initial states 



giving 
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A probability can be calculated, changing variables of integration such that 

^ ^ a and ^ /3, and evaluating with ^ = dn§ and ^ = dm^, 



Pi,r{e,N) = tip' 



,r+l 



(1 + eN) 



r+l+2 ■ 



The calculation for the intensity /(r) for p' proceeds as follows, 
/(r) = tr(a^ap") 

f , dO d<j) ( ( 1 \ ^] 
oc y dn— dm— |exp-(^e + ^j(n + m)| 



Now, 



= nil 
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+ fh\l + e'^\' + {..)e'^e-'^ + (..)e- 
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A. 2 Derivation of the visibility for thermal initial states 



and the latter two contributions integrate to 0. Hence 

I{t) oc |1 - e''^|^ydndmn''+^m'exp- l^e+i^ (n+m)| 

Evaluating the integrals is as for the probabihty calculation above. 



7(t) oc (|l-e'^|'(r + l) + |l + e^^|"(Z + l) 
oc I cos^ o + sin^ o + 



and the visibility is given by 



r + l + 2 
extremising at r = and t — tt. 
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B 

Derivation of the visibility for a Gaussian dis- 
tribution of the relative phase 



For a normal distribution with mean Aq and variance cr^, the probabihty density 
function takes the form: 



P(A) 



-(A-Ao)V(2a2) 



In what follows the visibility V is computed for a state pc of the form, 

-271 r^TT 
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27rP(A) + (contributions for 27r-periodicity in A) 

27r 27r 
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(B.l) 



where A = (J) — 6 and v iVe*^ denotes a Glauber coherent state with mean atom (or 
photon) number N and phase 9. It is assumed that the Gaussian width is small 
compared to 27r. 

The following identity can be used following the 27r-periodicity of the variables 
A, (f) and 9 and the lack of dependence on the mean phase variable 



2^ /-^^ d9 d<i) 
Jo 27r 27r 



^offset +7r 



^offset— 



27r 



(■) 



(B.2) 
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If a is less than 7r/2.58 ~ 1.2 then tipc = dA P (A) > 0.99 and pc is seen to 

be well normalised. 

V can be computed by comparing the probability for measurement at x, where 
the measurement is given by the Kraus operator M (x) = \/l/2N [e^'^^a + e~'"^^6), 
with the standard form 1 + cos (27rx — ip) (the momentum parameter k has be set 
to TT units) - refer Sec. 14.21 

This probability is computed as follows: 

tr [M{x)pgM{x)^] = /^""^rfA g-(A-A„)V(2<x2) ^ ^27rx - A)] 

= l + cos(27ra;-Ao) T"'"'-^ e'^^-^^^'/M cos (Aq- A) 
- sin (27rx - Aq) g-(A-Ao)V(2.^) (Aq- A) 

The second integral evaluates to because of change of sign about Aq. V may be 
identified with the coefficient of cos (27rx — Aq). Extending the range of integration 

to ±CXD, 

/OO -1 
rfA^^e-^'/(2'^')cosA 
-OO a\J1'K 

= e-^-^'. (B.3) 

This expression for V may be taken as exact when o" < 1. 

As an example consider the optical interference experiment analysed in Sec. 12.3.11 
of chapterElfor which two optical states, initially Poissonian with equal mean photon 
numbers, are combined at a beamsplitter. In the case that r photons are counted 
at the right output channel, the relative phase localisation is described by a scalar 
function c (6*, 0) oc cos^'^ ^ and the visibility V is calculated exactly as r/ (r + 1). 
For larger values of r c (6', 0) is well approximated by a Gaussian with standard 
deviation a = \/2/r and the calculation above suggests a value for V as exp (l/r). 
The fractional error in the latter estimate does indeed decrease for increasing r and 
is 0.26 for r = 1, 0.09 for r = 2, 0.04 for r = 3, less than 0.01 by r = 7 and less than 
0.001 for r = 23. 
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